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1 Abstract 

The effective residual interaction for a system of hadrons has a long tradi- 
tion in theoretical physics. It has been mostly addressed in terms of boson 
exchange models. The aim of this review is to describe approaches based 
on lattice field theory and numerical simulation. At the present time this 
subject matter is in an exploratory stage. A large array of problems waits 
to be tackled, so that known features of hadron-hadron interactions will 
eventually be understood in a model-independent way. The lattice formu- 
lation, being capable of dealing with the nonperturbative regime, describes 
strong-interaction physics from first principles, i.e. quantum chromody- 
namics (QCD). Although the physics of hadron-hadron interactions may be 
intrinsically complicated, the methods used in lattice simulations are simple: 
For the most part they are based on standard mass calculations. This chap- 
ter addresses commonly used techniques, within QCD3_,_]^ and also simpler 
lattice models, describes important results, and also gives some insight into 
numerical methods for multi-quark systems. 



*Prepared for World Scientific Publishing Company (International Review of Nuclear 
Physics, Vol. 9, Hadronic Physics from Lattice QCD, edited by A. M. Green) 
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2 Introductory Overview on Goals, Strategies, and 
Methods 



2.1 Modeling Nuclear Forces 

The theory of nuclear forces was pioneered by Yukawa j^] in 1935. A par- 
ticle of intermediate mass, a meson, accounts for the interaction energy of 
the proton and neutron. In the framework of a quantum field theory the 
meson appears as the quantum of an effective field that describes the strong 
interaction . The approximate range of the interaction 

fl=A (1) 

is roughly consistent with R k, 2 fm, for a meson mass in the region of 
m ~ 100 MeV. Those numbers are crude estimates based on the uncertainty 
relation (for energy and time) and the contemporary experimental evidence. 

An effective neutron-proton interaction can be obtained from solutions 
of the Klein-Gordon equation for the meson field (p 

{U + n?)ip{x) = gi;{x)ip{x) , (2) 

with the pion-nucleon coupling g. In the static limit, where the product of 
nucleon fields ^^ip is replaced with a generic point-like source, 

(_A + m2)(^(f) = 55(f), (3) 

the solution, i.e. the static Green function, 

-mr 

^{x) = j- (4) 

47r r 

is interpreted as a potential energy function, the 'Coulomb potential' of 
the meson theory. (Units are such that h = c = 1, conversion is he = 
0.197 GeVfm.) This line of thought is in close analogy to electrodynamics 
where the Coulomb potential emerges as the solution of a Poisson equation 
with a static point-like charge (source). 

Since the early work, the meson theory of nuclear forces has received 
much attention and refinements. One of the more significant insights [^J 
has been that the range of the nucleon-nucleon force can be subdivided into 
three regions: 

• The long-range region, inter-nucleon distances are above r « 1.5 fm. 
Here we have attraction. This is the region described by the Yukawa 
potential. 
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• Intermediate-range attraction, around r ~ 1.0 fm. 



• Short-range repulsion, below r ~ 0.5 fm. 
Other salient features are: 

• A tensor force, important for the quadrupole moment of the deuteron. 

• The presence of a spin-orbit force. 

• Charge independence of the strong interaction. 

In order to describe the short-distance region, scalar and vector mesons 
[r], p,uj, . . .) are needed in the boson-exchange model. Remarkably, a vr-vr 
correlated state (known as the a meson) is relied upon to produce interme- 
diate range attraction. The pseudoscalar meson (vr), which also gives rise to 
the long-range (though finite) Yukawa potential, contributes significantly to 
the tensor force. 

After approximately half a century the meson theory of nuclear forces is 
at a mature, sophisticated stage. Exhaustive review articles exist, for exam- 
ple Refs. @, ^, ^ 0] rnay serve the interested reader. Boson-exchange models 
of this type presently constitute the state-of-the-art for effective theories of 
nuclear forces. 

Typical for this class of models is an effective Lagrangian, say for the 
nucleon-nucleon interaction, such as 

s 



with pseudoscalar gp, scalar gs, vector gy and tensor /v couplings for the 
various bosons ip, and M is the nucleon mass 0. Potentials are obtained 
from tree-level diagrams to the nucleon-nucleon propagator. Those are most 
naturally written in momentum space. The simplest contribution derived 
from (^) is 



V{k) =- 



12M 



^ ^ ^ ^ ml Si2{k) 

(Tl • (T2 — 0"! • (72— + 



A;^ + nip k'^ + nip 



fi • T2 • • • , (6) 



where mp is the pseudoscalar mass. This already is the non-relativistic re- 
duction where the three terms correspond to a (coordinate-space) (5-function, 
a Yukawa, and a tensor potential, the latter containing 

Si2 = 3((Ti • k){a2 ■ k) - ai- (J2P ■ (7) 
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The momentum transfer is k = q' — q and ai, 02 are the nucleon spins. 
In addition, for each type of couphng momentum cutoff form factors are 
included with the couphng constants gp,gs, 9v and /v , for example 



9l 



9l 



A2 



mi 



(8) 



The inverse of a momentum cutoff parameter Ap is a length, and as such 
indicates the region of validity of the model. Values of ~ 1 GeV, or equiv- 
alently ~ 0.2 fm are typical. Other refinements of the model include isobar 
(A) degrees of freedom, and the extension to the strange hadron sector j^, ^. 
Typical potentials are illustrated in Fig. |^. 

300 I 
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Figure 1: Generic form of a central nucleon-nucleon potential according to 
boson exchange models and its contribution from various bosons. Schematic 
drawing after Ref. 0]. 

The philosophy of boson exchange models is to take the phenomenology 
of low-energy excitations of the hadronic vacuum as a starting point. The ex- 
perimentally known spectrum of mesons and hadrons is used to reconstruct 
the forces between those particles. The result is an effective theory with ad- 
justable parameters and predictive power in a reasonably well-defined realm 
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of validity. An example is shown in Fig. ^. Naturally, the question of how 
hadronic interaction results from the underlying microscopic structure of 
hadrons and the, as we now know, complicated structure of the hadronic 
vacuum, is beyond the scope, and aim, of the model. 




0.0 0.4 0.8 1.2 1.6 2.0 4 8 1.2 1.6 2.0 

P,,„ [GeVVc] P,^^ GeV/c 



Figure 2: Nucleon-nucleon scattering phase shift for channels ^Sq and ^Z?2 
compiled from Ref. [10|. In the channel the smooth curve is a 'Bonn- 
potential' result. 



In the present context we wish to look at hadronic interactions from 
the perspective of lattice QCD. Since the physics of hadronic interactions 
is contained in the low-energy, non-perturbative, regime of QCD the lattice 
formulation is well suited for this task. 



2.2 The Lattice QCD Perspective 

From the lattice QCD point of view the three-fold subdivision of the inter- 
nucleon interaction region, see Fig. ||, appears in a new light: 

• The long-range region is expected to be dominated by excitations of 
the QCD vacuum that are recovered only in the chiral limit of the 
lattice simulation. If Wilson fermions are used this region can be 
monitored by extrapolating to the zero-mass limit of the pseudoscalar 
meson, which is the Goldstone boson of QCD. The mechanism for pion 
exchange is available on the lattice because correlated quark-antiquark 
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pairs are dynamically provided by the lattice action. This is true even 
if the simulation is limited to the quenched approximation. 

• In the intermediate-distance region various excitations of the QCD 
vacuum of somewhat higher mass will take over. Those may be pure 
hadronic states, gluonic excitations, or combinations of such. In con- 
trast to boson exchange models the precise nature of those need not 
be known as input. Again those would be accounted for by the lat- 
tice action. At r ~ 1.0 fm this is really a transition region between 
an appropriate description in terms of effective hadronic fields and 
the asymptotically free regime of the interior quark-gluon cores of the 
hadrons. The validity of traditional boson exchange models seems 
somewhat stretched here. With the currently available lattice con- 
stants (~ 0.1 - 0.4 fm) the intermediate-distance region appears ap- 
propriate for lattice QCD. 

• For small distances the compositeness of the hadrons, as finite-size 
quark-gluon systems, will play the dominant role. Not surprisingly, 
quark models have been used to study the nucleon-nucleon interac- 
tion. The energies involved are naturally larger, for example « 2 GeV, 
and thus point to lattice constants of about ~ 0.1 fm, being not out 
of reach of contemporary lattice simulations. Lattice QCD provides 
dynamically for the compositeness of the hadrons. Internal hadronic 
structure is, in turn, crucial for the features of the residual interaction. 

Thus lattice QCD is in a position to combine the classic nucleon-nucleon in- 
teraction regions in a single unified approach. Moreover, the strange hadron 
sector, or effects of strangeness in the light-quark sector are straightfor- 
wardly incorporated. 

2.3 Short and Long Term Goals for Lattice QCD 

While boson exchange models have a long tradition, attempts to study resid- 
ual interactions between systems of hadrons within lattice QCD are only in 
an experimental stage. The physics of hadronic interactions in terms of first 
principles is rather involved, see Fig. ^. Conceptual, technical, and com- 
putational questions are very much in the foreground at this time. The 
tentative nature of this situation will be apparent throughout this chapter. 

Nevertheless, there is potential for new physics, or at least for under- 
standing known phenomena in a new light. In particular, probing hadronic 
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Figure 3: Cartoon of a strongly interacting hadron-hadron system from the 
QCD perspective. 

matter at somewhat higher energies, say in the 4-8 GeV region with elec- 
tron beams, has become a new focus in experimental nuclear physics. The 
interest here lies in the boundary between particle (high energy) and nu- 
clear (low energy) physics which we will refer to as hadronic physics. Our 
short-term working list is as follows: 

• Discuss various basic concepts and numerical methods using simple 
lattice models. 

• Study simple systems, mostly meson-meson, and extract non-relativ- 
istic potentials. 

• Explore alternative approaches to the problem, like momentum space 
and coordinate space formulations. 

As a long-term perspective we have in mind: 

• Study realistic meson-meson, meson-baryon, and baryon-baryon sys- 
tems. For example, the pion-nucleon system features the prominent 
delta resonance, a hallmark of strong-interaction physics. 

• Residual interactions in the strange-hadron sector. 



• Develop advanced methods, possibly the direct extraction of scattering 
amplitudes from the lattice. 



2.4 Probing the Lattice 

A hadron-hadron force may be viewed as a hyperfine interaction in the 
hadron-hadron system. In one way or another, we must therefore compare 
the mass spectrum of the interacting hadron-hadron system to the spectrum 
of the non-interacting hadron-hadron system. 

About two decades ago it was reahzed that finite volume effects may 



in principle be exploited to achieve this end |11]. The energy levels of a 
two-hadron system, say two mesons with mass m each, are subject to finite- 
size corrections. Those decrease with powers of 1/L as the lattice size L 
increases. The formal study of finite-size effects in this context is due to 
Liischer |l^, 13 1. In Ref. [|^] the following equation for the ground-state 



two-body energy shift 

Wo-2m = {l + + + '^(^"'^ 

is derived. It is valid for s-wave elastic scattering, below the particle pro- 
duction threshold. Above ci = —2.837297, C2 = 6.375183, and oq is the 
s-wave scattering length. This result has been the basis for a number of 
studies (H, H, m, 0, H, 1^, P, H, m, H, H, m. The extraction of 



scattering information beyond the static limit, like phase shifts at nonzero 
relative energy, is formally more involved [^, |2^ . 

In general it will be necessary to probe the lattice (vacuum) with one- 
and two-hadron operators. Let (j)j;{t) be a composite field, say, constructed 
from quark-antiquark (meson) or from three-quark (baryon) fields. Let p 
be the hadron momentum, for simplicity we disregard additional quantum 
numbers. The time correlation matrix built from one-hadron operators 

4?(*'*o) = {^{t) c^^ito)) - {4{t)){c^^^{to)) (10) 

(2) 

possesses eigenvalues Xrn{{t-,to) which behave exponentiahy, 

A(^) (t, to) oce-"""! (11) 

The time correlation matrix is just the euclidean version of a propagator, 
in our example a 2-point function C^^^ on the hadronic level (a 4- or 6- 
point function on the quark level). Thus from ( pT] ) the energies Wmi of the 
one-meson system can be extracted. Typically, this is done numerically. A 
practical example is shown in Fig. ^. The interpretation of the is that 
those are energies of the hadron moving through the lattice. In principle, 
the hadron could be in an internally excited state. 
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Figure 4: Typical time correlation functions (eigenvalues) extracted from a 
lattice simulation. The examples are for a pseudoscalar meson (PSM) and 
a vector meson (VEC) using QCDg^^ and Wilson fermions. 

The fields ^p{t) = (t>-p{t) (t>+p{t) can be used for probing two-hadron 
states. These are operators with total momentum P = 0. Note that p in 
this example means the relative momentum of the two-hadron system. We 
construct a 4-point correlation matrix 

Cflit, to) = $,-(*o)) - (4(0)(^,-(io)) , (12) 

and, correspondingly, the eigenvalues 

A(,4)(t,io)oce-^"(*-*°) (13) 

yield the energy levels Wn of the two-hadron system. Most importantly, 
these levels somehow contain the desired residual interaction. 

A sensible definition of the latter, of course, also involves the eigenvectors 
of C^^) and C^^\ however, comparison of the level schemes 

m={mi,m2) = Wrm+Wm2} to {Wn} (14) 
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already gives important insight into the nature of the hadron-hadron force. 

Probing the lattice with momentum-type operators is only one of several 
alternatives. The use of coordinate space operators is another popular choice 



3g, H, H, |4|,|3|. 



2.5 Finite-size Methods 

The spectrum {VFn} of the interacting two-body system may be utilized in 
a variety of ways. Liischer has shown how to extract scattering phase shifts 
directly from the spectrum The idea is to match the poles of (inter- 
acting) Green functions that live on the finite-sized periodic lattice to the 
poles of (free) Green functions in infinite volume. Since the lattice Green 
functions comprise the interaction, comparison in the asymptotic (free) re- 
gion yields the phase shifts. Using methods similar to the 'Jost' formalism 
of scattering theory [ 36 1 an equation of the following form 



det 



2i5 



{M + i)-\M -i) =0 (15) 



is derived. Here M = M{q^) is an analytically constructed energy-dependent 
matrix iW"^ = m?+q^). It reflects features of the 'natural' (free) spectrum of 
the periodic lattice. The other matrix e^^^ in (|l5|) contains a set of scattering 
phase shifts 5i{q^). The lattice simulation gives a discrete spectrum {VF„}. 
A certain energy Wn, for a fixed n, is then used to calculate M{q^) and ( p^ ) 
is solved for the corresponding 5i{q^). This yields a set of phase shifts, as 
many as the size of M, at each energy Wn- The physical 5^{qf^), at a discrete 
q^, is contained in this set. 

A simple illustration of the idea behind this formalism is the follow- 
ing: Consider a solution ip{x) of a Schrodinger equation in the noninteracting 
case with a phase kL = 2tt n. It satisfies ip{x) = ip{x+L), see Fig. |5|(a). If the 
interaction is turned on, say an attractive potential with range |x| < 0.2L, 
the periodic boundary condition would be violated, see Fig. |5|(b), unless 
compensated for by a simultaneous change in the kinematical phase kL. 
In other words, the change in kL necessary to restore the original periodic 
boundary conditions, see Fig. ^(c), is a measure for the phase shift. 

The kinematical phase kL is changed by varying k ox L. Varying L 
is reminiscent of a finite-size effect, see (|9|). Finding the change in k at 
fixed L is equivalent to finding all (discrete) two-body levels for the 

interacting case. This method, which already appears in Ref. [O], has been 



used by Liischer and Wolff |27] to obtain scattering phase shifts in an 0(3) 



nonlinear sigma model on a 1-1-1 dimensional lattice. 
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Figure 5: Generic solutions ^l^{x) of a Schrodinger equation: (a) no interac- 
tion, (b) an attractive potential around x 0, and (c) with decreased wave 
number k to restore the original phase conditions. 

2.6 Residual Interaction 

Scattering phase shifts are very close to observable quantities (cross sec- 
tions). This is an appealing feature of Liischer's method, as the above 
procedure came to be known. On the other hand an effective residual inter- 
action between hadrons is more in line with the historic development and 
has a heuristic value which should not be underestimated. We will therefore 
place considerable emphasis on this aspect. In this approach the informa- 
tion contained in the 2- and 4-point correlation matrices, see ( [lO| ) and ([T^, 
will be utilized in a different way. 

In principle an effective residual interaction is subject to definition. To- 
wards this end we will first split a free, or non-interacting, part off the 
4-point correlation matrix (|l^) 

C<^\t,to) = C{t,to) + Ci{t,to), (16) 

where C is built from 2-point correlators C^^^ and describes two non-interacting 
hadrons on the lattice. This can be accomplished by means of diagrammatic 
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considerations. In terms of (16) we then define an effective residual interac- 
tion as 

In principle the interaction Hamiltonian (potential) is non-local. Our focus 
is to compute Hi from a lattice simulation and eventually calculate the 
corresponding phase shifts. 

The above method relies on lattice sizes large enough to accommodate 
two hadrons. Assuming a manageable lattice x T with, say, L = 10 
and physical extent of about 4.0 fm this translates into a lattice constant 
of a = 0.4 fm. As experience with hadron spectroscopy has shown the 
standard Wilson plaquette action requires a ~ 0.1 — 0.2 fm in order 
to work reasonably close to the continuum limit. This does not sound very 
encouraging for putting multi-hadron systems on the lattice. 

2.7 Use of Improved Actions 

Fortunately, considerable progress in recent years has provided us with so- 



called 'improved' lattice actions |38|. Here, in the simplest form, one supple- 
ments the Wilson action, built from elementary 1x1 plaquettes (pi), with 
larger 1x2 rectangles (rt) 

Sg[U] = 13 {Spi[U] + CrtSrtP]) . (18) 

The Wilson part of the action (pi) approaches the classical (tree-level) limit 
like o(a^) as a ^ 0. If Crt = —1/20 the (rt) term in (|l8|) improves this be- 
havior to order o(a^). We refer to such an action as o(a) tree-level improved. 
Better yet, one can also account for a large part of the quantum fluctua- 
tions of the gluon field through renormalizing the link variables U^{x). Of 
course, this would destroy the tree-level cancellations between the (pi) and 
(rt) terms in (|^). Within mean-field theory it can be 'repaired' by changing 
the factor Crt to 

^'■* = -2^' 

where u, known as the tadpole factor, is the fourth root of the average 
(4-link) plaquette 

u = {hieT,Upi)^'\ (20) 

The tadpole factor can easily be computed self-consistently within a run- 
ning lattice simulation. There are also radiative corrections which can be 



incorporated into the lattice action to some extent [39, 4C, 41 1. 
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The benefit of improvement is that the physical continuum hmit, a — > 0, 
is approached 'faster' compared to the Wilson action. The 'same physics' is 
described with improved actions at larger lattice constants a. For example, 
a Wilson action at o ~ 0.1 fm may give a similar mass spectrum as an 
improved action at o 0.3 fm. For our present purpose of putting two 
hadrons on the lattice, this means that, stretching things a bit and choosing 
a ~ 0.4 fm, we would only need an L'^ x T lattice with about L = 10 for a 
lattice of 4.0 fm across. This should give us a first window on two-hadron 
systems on the lattice. 

Improved actions have been a big innovation for lattice field theory |42]. 
There are versions of improved actions including fermions and work 

pushing the improvement to higher order, for example o{a?) improved ac- 
tions |44, 46|. - In Appendix |8| we have collected some topics that 
are relevant to this chapter. - Strong-interaction physics, in particular the 
boundary between particle and nuclear physics, is now open to studies within 
quantum chromodynamics. 



3 A Simple U(l) Lattice Model in 2+1 Dimen- 
sions 

In order to discuss some practical work we look at a simple but non-trivial 
model. We consider U(l) gauge theory on an L'^ X T lattice, in d + 1 
dimensions, with d = 2. The intent is to study the interaction between 
two mesons living in a plane. For the sake of simplicity we work with the 
(unimproved) Wilson action and use staggered fermions. 

3.1 Lattice Action 

For a U(l) gauge group the link variables are 

U^{x) = e*®'^^^) G U(l) , (21) 

with fi = l...d+l, X = {x,t), and Q^{x) £ [— 7r,7r). The plaquette (pi) 
variables are 

Upi{fi,i^;x) = U^{x)Ui,{x + fi)Ul^{x + i')Ul{x) 

= exp (iQpi{^,u-x)) , (22) 

where 

Qpiin, V] x) = Q^{x) + 9^(2; + fi) - e^(x + v) - e^(x) (23) 
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is the oriented plaquette angle. We write jl for a vector of length a in 
direction ^u. The (compact form of the) Wilson gauge field action is 



X ll<V X 

(24) 

For simplicity, fermions are treated in the staggered, or Kogut-Susskind, 



scheme where Dirac indices have been 'spin-diagonalized' away |47, 48, [l9f| . 
Thus we work with one-component Grassmann fields Xf{x) defined on the 
lattice sites x = (x, t), with (external) flavor index f = u,d, s. There is just 
one color index, C = 1, which we omit. The fermionic action can be written 
as 

^^y f 

with the flavor-independent fermion matrix 

G'^ix,y) = ^J2'^t'^^) Ut^i^)^x+fi,y - Ul{y)6x^y+f, +mFSx,y, (26) 

including a mass term. Here r]fj_{x) = (_i)^i+---^m-i , ^ = 2 . . . d+1, r/i(x) = 1 
is the staggered phase. 

This lattice field model is confining in finite volume [50], and exhibits 
chiral symmetry breaking [^] and monopole condensation [p^ ], for appropri- 
ately chosen inverse gauge coupling /3. Its basis is the euclidean generating 
functional integral 

Z[J,J]= j[dU][dx\[dxK^''^^^'^^^^''>''^''+^''^^\ (27) 

with the fermionic sources Jx = Ylix J{x)x{x) and xJ = x{x)J{x). We 
will use gauge configurations in the quenched approximation. 



3.2 Meson Fields 

The simplest hadrons we can work with are pseudoscalar mesons. In the 
staggered scheme suitable one-meson operators are 

<i>^{t) = L-2 e'^-'^Xdix, t) Xu{x, t) . (28) 

X 

The X sum extends over the spatial sites of the lattice. The fixed flavor 
assignment u,d leading to a vr^ will simplify computations. For d = 2, a 
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planar lattice, the momentum parameters are 

p=^{ki,k2), where A:i,2 = -(^ - 1) • • • ^ , even L , (29) 



if L is odd we have ki = -{L - l)/2 . . . (L - l)/2. The operators (|2|) 
probe the lattice (vacuum) for composite meson states with momentum p. 
Products of those operators probe for multi-meson excitations. For example 

^p{t) = (t^^p{t) Mt) (30) 

describes two-meson fields with total momentum P = and relative mo- 
mentum p. The reduced mass is m/2. 

3.3 Correlation Matrices 

The time-correlation matrix for the one-meson system is 

4?(*'^o) = (4(i)0,-(to)) - (4(t))(0,-.(to)>. (31) 

It can be worked out with Wick's theorem in terms of contractions between 
the Grassmann fields. In our simple example we use, see Appendix |^ 

...Xf {x)xf' {x')... = ...5ff.G{x,x')... (32) 
...%{x)Xf> {x')... = ...6ff,G*{x,x')... , (33) 

where n indicates the partners of the contraction, and G is the inverse 
fermion matrix, cf. (^). This may be used to work out the correlation 
matrices. For example, consider 

'^^-(*o)> = E E e-'^^-^e^'^-HxuiS, t)xd{x, t) Uv, to)Xu{y, to)) • 

X y 

(34) 

Because of the flavor assignment in (^) the separable term in ( |3ll) is zero, 
and there is only one group of contractions in the first term 



with 



2 1 1 2 , ^ 

C^^' ~ XuXdXd^u. (35) 



1 1 

■■■^d{x,t) Xd{y,to)--- = ■■■G{xt,yto) ... (36) 
... |„ (x,t) X„ (y,to)--- = ...G*{xt,yto)... . (37) 
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Thus we obtain 

c(5(t,to) = L-^X^X]e-^-+^'^-^l|G(xi,yto)|2) . (38) 

X y 

Assuming translational invariance 

(|G(f,t;y,to)|') = {\G{x + a,t;y + a,to)\^) (39) 
renders the 2-point correlator diagonal 

4?(*'*o) = '^P?^"'E^"''"^l^(^*'^o*o)|') e*^^-^^°. (40) 

X 

The point xq is arbitrary, C^"^^ is independent of xq, and the phase factor is 
irrelevant. 

The 4-point correlator describes the propagation of two interacting me- 
sons on the lattice 

C^i(i,to) = (4(0 Mto)) - (4(t))(<I>,^(to)) . (41) 

Here p and q are the relative momenta in the meson-meson system. In terms 
of the fermion propagator this correlator reads 

ci^l{t,to) = ^^Y^Yl e^P-(^2-xi)+ig-(y2-j/i) (42) 

xi X2 yi y2 

{{G*{x2t,y2to) G{x2t,y2to) G*{xit,yito) G{xit,yito) 

+ G*{xit,y2to) G{xit,y2to) G*{x2t,yito) G{x2t,yito) 

- G*{x2t,yito) G{x2t,y2to) G*{xit,y2to) G{xit,yito) 

- G{x2t,yito) G*{x2t,y2to) G{xit,y2to) G*{xit,yito))) . 

For an SU(N) gauge group, G Gcc'-, there would be sums over four sets 
of color indices corresponding to the xi,X2,yi,y2 assignment in (p2|). A 
diagrammatic classification of the various terms which contribute to C*-^-* 
proves useful. Such a classification arises naturally from working out the 
contractions between the quark fields in (pl|) and (31). Let us write 

(27(4) ^ fjiAA) _^ fj{4B) _ (j{4C) _ (j{AD) (43) 

+ X - M - M (4^) 
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for the four terms as they occur in (E|). Each dia gram hne corresponds to 
a fermion propagator, see Fig. Specifically, using the notation introduced 
m (H) and (ID, we have 



43 21 12 34 



43 34 21 12 

■At , 



(45) 



21 43 12 34 



43 34 21 12 



(46) 



where the n = 1 ... 4 identify the partners x and x of a contraction. In 
diagrams A and B no quark exchange between the mesons occurs. Diagram 
B describes the exchange of (composite) mesons as a whole. On the other 
hand diagrams C and D both contain quark exchange between the mesons 



23 41 12 34 



(47) 



41 23 12 34 

= {cl>l^^l^^_^^^.) . (48) 
Thus these diagrams must be considered as sources of effective interaction. 

(xt) (xt) (xt) (xt) (xt) (x,t) (xt) (xt) 



(y,g (y,g (y,t„) 






(y,g 



A 



B 



C 



D 



Figure 6: Diagrammatic classification of the correlation matrix C^^^ accord- 
ing to (H). 

The propagation of free, noninteracting, mesons is contained in C^'^^^ + 
However, gluonic correlations between the mesons are also present 
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because the gauge configuration average (...) is taken over the product of 
all four fields. (This is also true for C^^^^ and C^^^\) Gluonic correlations 
of course contribute to the effective residual interaction. Hence it is essential 
to isolate the uncorrelated part contained in C^^^^^ + C^^^^ . 



3.4 Correlation Matrix for Noninteracting Mesons 

Let us examine C^"^"^-* . Expressing the contractions in ( ^5[ ) in terms of quark 
propagators with the notation introduced in ( |32[ ) and (^) we have, in generic 
notation 

/ 4 3 2 1 

~ (G*GG*G), (49) 



where the ~ indicates the Fourier sums, etc, which carry over from (p8[), see 
(^2|). The gauge configuration average in ( ^9|) may be analyzed by means of 
the cumulant (or cluster) expansion theorem^. Taking advantage of {G) = 
we have 

4 3 2 1 

{G*G G*G ) = (50) 

43 21 41 23 42 31 4321 

{G*G ){G*G ) + {G*G ){G*G ) + {G*G*){G G ) + {{G*G G*G )) . 

The last term defines the cumulant. The first three terms on the right-hand 
side of ( pO[ ) are illustrated in Fig. |^ The dashed lines indicate that the 
propagators are correlated through gluons. Evidently only the first one of 
the three separable terms in ( pO[ ) represents free, uncorrelated, mesons. All 
other terms in (^0|) are sources of residual effective interaction between the 
mesons. We therefore define 



4 3 2 1 



43 34 21 12 

C^'""^ = {^lr/P+a)i^-J-a) ~ ) {G* G ) . (51) 



A similar analysis of C'-'^^-* leads to 

43 34 21 12 



The sum of those is the free meson-meson time correlation matrix 

C%{t, to) = C^^^^ {t, to) + C^^P {t,to). (53) 



^The cluster expansion theorem involves the definition of the cumulant as a general- 
ization of the standard deviation. For two random variables Xi and X2, for example, we 
have {X1X2) = (Xi)(X2) + ((X1X2)). If Xi — X2 the square of the standard deviation is 
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<G'G><G'G> = 





<G*G*><GG> = 




Figure 7: Illustration of the cumulant expansion of C^^^ . The dashed lines 
represent gluonic correlations. 

— (4) 

C is an additive part of the full 4-point correlator, 

CW = c(^)+C7f). (54) 

The remainder 

Cf = + - C^'^) - - (55) 

contains all sources of the residual interaction, be it from gluonic corre- 
lations, or quark exchange (or quark- antiquark loops, if the simulation is 

unquenched) or any other effects. The free correlator C^^'' describes two 
noninteracting identical lattice mesons. Their internal structure and masses 
consistently arise from the dynamics determined by the lattice field model 
and its numerical implementation. As already mentioned the separable term 
in ( |3T[) is identical zero because of the u, d flavor assignment to the Grass- 
mann fields of the one-meson operators (p8|). 

C^^^ may be expressed in terms of the 2-point correlator C^^\ Using 
(|l|) and (|5l|),(|5|) gives 

-^(4) _ (2) ^(2) ^(2) ^(2) , . 
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C^^^ is diagonal in the momentum indices. Writing 

where Cp'(t, to) is given by ( |40| ) and has the obvious property 

c_j?(t,to) = 4(i,to), (58) 

we obtain 

Cp]{t, to) = + 5-p-,,-) I cp-(t , to) I ' . (59) 

To be consistent with the computation of cj^l{t,to), using (|4^), Cp'(t, to) 
should be computed from (^). It is interesting to note that the property 

77(4) _ 77(4) _ 7^(4) _ 7=7(4) 

which is evident from (|56| ) and is explicit in (^2|), also holds for C^^^. This 
property reflects Bose symmetry with respect to the composite mesons. Per- 
mutation of the mesons, one with momentum +p the other with momentum 
—p, results in the substitution p —p. Clearly, combining diagrams C^^"^^ 
and C^^^^ into C^^^ is crucial for the symmetry ( |60|) to hold. 

3.5 Computation with Random Sources 

The computation of C^^^ from the expression (^) is a formidable numeri- 
cal task. First, propagator matrix elements are needed for arbitrary y, as 
opposed to only one column of G like in (^0|). Typically G is computed by 
solving a linear equation like G~^X^^^ = R where i? is a given source vector; 
the solution vector X^^'^ then represents 'some column' of G. For arbitrary 
y such a strategy translates into L*^ linear equations, all for a fixed time slice 
to and not even counting color indices (or Dirac indices in case of Wilson 
fermions) . 

Furthermore, the sum over lattice sites in ( ^2|) contains L^'^ terms. Even 
for d = 2 and our simple U(l) gauge model those are ~ 10^ terms. In a 
realistic calculation, for d = 3 and SU(3), these would be L^'^S^ 10^^ for 
a meson-meson system, and even larger if baryons are involved. It is clear 
that an efficient numerical technique is called for to deal with computing 
quark propagators and correlator matrices. 



A practical way is the use of random estimators [54, M, 56, p^]. Let 



R{x) be a vector of independent complex random deviates with 

{R{x)R*{y))R = S^,y, (61) 
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where (. . denotes the random R-average. Gaussian deviates come to 
mind, but Z'^ noise also works and may sometimes be more desirable |58]. 
The idea is to solve 

Y,G-\x,y)X^^\y) = R{x) (62) 

V 

for ^ recalling that the fermion matrix G^^ is known, then the statistical 
average 

^X^''\y)R\z))n = G{y,z) (63) 

is an estimator for the inverse G of G^^ . 

Since random sources K with various characteristics may be chosen the 
technique is very versatile. Let us expand our notation, where 

[x) , with r = 1 . . . AT^ , (64) 

is a finite sample of vectors drawn from the distribution. The label a 
distinguishes between various types of sources, for example we take a = 
to mean 

Here, random sources are placed on one time slice only. A more complicated 
example is a = {to,p) with 

^(r;to,p)(^^ t) = e'P-^R'-'''^°'P^{x)6l° . (66) 

These are Fourier-modified (boosted) sources on time slice to- 

Given a finite sample (|64|) , an estimator for the statistical average of (|6l|) 

is 

_±_ Rir;»)(^x^ji{r;f^r^y-j ^ i?(^;°)(x)/?('^''3)*(2/) = . (67) 

^ r=l (r) 

The above, intuitive, notation ^^r> ^'-'^ statistical R-average (. . is 
convenient in expressions of correlation functions. The second 5 on the 
r.h.s. in ( p^ ) comes from drawing stochastically independent samples for 
each type of source. 

For definiteness we continue with sources of type ( |6^) set at a fixed time 
slice to- Then the solutions X of the linear equations (|62|) 

J2 G-\xx,,yy,)X(^-''o){yy,) = i?(^;*°)(f)4» (68) 
yy4, 
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can be used to estimate columns of the propagator G. Applying G to both 
sides of (^) gives 

X(^'*°)(zz4) = Y,G{zz4,xto)R^'^''''\x) . (69) 

X 

Multiplication with R^'^'^^'^* (y), taking the R-average and using (^) yields 
G(zz4,yto) = j;x('-;*«)(zz4)i?(^'*°)*(y) . (70) 

(r) 

In this form the random estimator can be used directly in the expression 
(ID and (g|) for C(2) and C^^), respectively 

Another, potentially useful, relation follows from (|69|) 

^X(^^*°)(zz4)x('-*)*(yy4) = G{zz,,xto) G{yy^,xtor . (71) 

(r) X 

A more general version of ( [7l| ) can be derived when Fourier-modified sources 
(|66|) are used. We have 

Y^X^^-''^'P^{zz^)X^^-''^^'^*{yy,) = L-'Y.^^^^ 
(r) X 

(72) 

Direct application of ( [72D gives the 2-point correlator ( |38D in the form 

C^J{t, to) = {L'^ e-'^'^Y X^"'''"'^{xt) X('^'*0'°)*(xt)) . (73) 

X (r) 



This equation is particularly useful for testing diagonality, see (4C). 

In the case of the 4-point correlator two independent sets of random 
sources are necessary. The four contributions C^'^^^ through C^'^^\ see (^3|), 
may be computed separately, for example 

C^/\t, to) = (YYl X^^'''°^'^\xit) X('"i;*«'^^*(i;it) 

(n) XI 

(?-2> X2 

There are similar expressions for and C^^^), C^^^). The noninteract- 

ing correlator C^^^^ is computed via (|73|) and the use of (|56|). 

Having selected an even number A'^/j of random sources, in the above 
case a subset of Nr/I of those can be used in Yl{ri) other half in 

Yl{r2)' There are (i^) possibilities to choose such a subdivision (e.g. 70 
for Nji = 8 and 12870 for Nr = 16). 
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4 Effective Residual Interaction 



Potentials are not unique in the sense that a class of potentials can produce 
the same phase shifts (phase-equivalent). Therefore an effective residual 
interaction extracted from the lattice is subject to definition. In this section 
we study the perturbative expansion of the time correlation matrix for an 
interacting elementary Bose field on a euclidean lattice. Our intention is to 
find guidance towards a sensible definition for the lattice simulation. 

4.1 Perturbative Definition 

Let Co = Co{(j),d(l)) be the free Lagrangian for an elementary Bose field 
(j){x) defined on the sites x = {x, t) of the lattice. It is understood that (p 
is subject to the usual canonical quantization through commutators. Let 
£ = £0 + Ci such that Cj = Ci{(j)) is a (small) interaction. In the usual way 
£ gives rise to a Hamiltonian 

H = Ho + Hi, (75) 

where Hq is the free part and Hj a perturbative interaction. In view of ( |28| ) 
and (l30|) define 

X 

and 

^p{t) = 4>-p{t) hpit) ■ (77) 

In the correlation matrix 

cfi{t,t,) = ml{t)^<i{tcm (78) 

the separable term of ( pl) h as been dropped since it is zero for the fiavored 
quark fields, like in Sec. |3.3| . The (nondegenerate) vacuum state |0) satisfies 
H\Q) = Wo|0). We will assume that its energy is zero, Wq = 0. Thus the 
time dependence of the correlator ([T^) may be made explicit 

cf.{t,Q = (0|dt.(t,)e-^(*-*^)c|,KO|0). (79) 

Switching to the interaction picture, we define 

Hj{t) = Hi e"^" (80) 
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and the (euclidean) time evolution operator^ 

U{t, Q = e^«(*-*^) e-^(*-*=) . (81) 
The perturbative expansion of the latter 

U{t, Q = / dti... dtN T[Hi{ti) . . . Hi{tN)] (82) 

then induces a perturbative expansion of the correlator 

C^%Q = (0|<i;.(Oe-^o(*-*^)C/(t,te)l>,-(tc)|0) (83) 

oo 

= E^"^(Mc). (84) 

N=0 

The zero-order term (7(4;A^=o) evidently describes noninteracting mesons. 
Order N = 

Let \w) be a complete orthogonal set of eigenstates of Hq 

Ho\niy) = wjP^\niy) , (85) 

where wji^^ are free two-meson energies on the lattice, and is a degeneracy 
index. Lattice effects set aside, those energies should be close to m? + 
with m being the rest mass of one meson and p = \'p \ its momentum. Also, 
define the relative meson- meson momentum-space wave functions ipnu{p) 
through 

4'lVi°J(p) = (H^'p(*c)|0)% (86) 
where Cnv are normalization factors. The order = correlator then is 

cl^='\t,Q = |cW|^e-^'^°'(*-*=)^W(p)^(t(?) • (87) 



With properly chosen normalization factors Cnu we expect orthonormality 
and completeness 



This is in analogy to the transfer matrix formahsm. 
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pq ■ 



(89) 



For a free elementary Bose field this is almost a trivial point since the ipnu 
are merely plane (lattice) waves. A glance at ( |87|) shows that, technically, 
those could be obtained as (normalized) eigenvectors from diagonalizing the 
correlation matrix (7(4;A^=o)^ where Icn'jyp are the eigenvalues, at t = tc- For 
the case considered here all eigenvalues of (7(4;^=o) -^^iu \)q nonzero. 

In a lattice model, however, where the role of $p'(tc) is assumed to be 
a composite operator made from fermion fields, see (p^) and (|30|), it can 
not be a priori excluded that an operator matrix element, of the type as it 
occurs in (p6|), is identically zero for all p. In this case the free correlator 

C^^^ would have an eigenvalue zero (for all t). Likewise, if, for some reason, 
the set of hadron operators used to construct the correlation matrices on the 

— (4) 

lattice is linearly dependent, one should expect C to have an eigenvalue 
zero. 



Order = 1 



From (p^)-(|8^) we obtain 



pq 



-Ho{t~t.) fdhHi{h)%iQ\0). (90) 

Jtc 



Upon inserting the complete set \niy) on both sides of Hj{ti) and using (|80| ) 
the ti integral over exponentials can be carried out explicitly. The result is 



-EE*'™*?)*™",''*'?) 

(nz.|^,|m;u)cW*cW exp 



sinh 



(t tc)5niji + 



2 

{t - tc) 



(1 - 5„ 



(91) 



Without loss of generality the normalization constants Cnu may be chosen 
real and positive, with the phase factors being absorbed into ipnu{p), as is 
evident from ( p6| ) . A glance at ( p7| ) shows that the two normalization factors 
and the exponential in ( |9T|) may be removed by multiplying the correlation 
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matrix c(^'^ ^) from both sides with the inverse square root of C^'^'^ ^\ 
Hence the matrix elements of 

^(4;^=i)(t,t,) = (7(4;^=o)(t,t,)"'/'(7(4^^=i)(t,i,)c(4^^=o)(t,f,)-'/' (92) 

in the basis ipnuip) are products of {nulHjlmfi) and the expression inside 
{. . .} of (|9l|). The t derivative of the latter is equal to one at t = tc- Thus 
we have 

= - E E V'i°J(p)(nHi?/|m/.)VLt (9) ■ (93) 

t—tc 

Using (|8^ ) and (|89|) , this translates into an explicit equation for Hj inde- 
pendent of the basis 

(94) 

which is valid to order = 1. Finally, we may replace (7(4;^=i) in (|9|) 
with the full correlation matrix C^^\ The corresponding expression ( p^ ) 
for Hj will still be valid up to order = 1 in perturbation theory. Thus, 
summarizing, define the effective correlator 

^W(t,0 = C'(^^^=o)(t,t,)-'/'cW(t,OC(^^^=°)(t,te)"'^' (95) 
understood as a matrix product, then the meson-meson interaction Hj sat- 

IsflGS 

e:W(t,t,) = e-^^ (96) 

up to (at least) order = 1 in perturbation theory. 

The utility of these results in the framework of a lattice simulation lies 
in the analogy which can be drawn between (7{4;^=o) ^-^^^ ^j^g £j.gg correlator 

C^^\ and between (7^^^ and the full correlator C^'^\ The analogue of (]96| ) 
may then be considered as the definition of an effective interaction. 

Order = 2 

A calculation similar to the above reveals that (^) is valid to order A^ = 
2. However, for N > 2 the situation is complicated by the presence of 
disconnected diagrams, see Fig. |8[ The diagramatic series of connected 



dt 



Hi 



dt 
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insertions of Hj still gives rise to the exponential form 
that, within the region of validity of (|96|), the operator 



Hi 



51nCW(t,fc 
Ft 



is independent of t and tc 



1). We conclude 



(97) 



C) 



o 



o 

i) 



Figure 8: Low-order diagrams for the elementary meson field correlator. 
The circles represent Hj. 



4.2 Effective Interaction for Composite Operators 

The above result is now easily transcribed to a system of composite hadrons. 
In this case the effective correlator (^) becomes 

= C^''\t,tof'^' ci''\t,to)C^''\t,tof'^\l (98) 

4'^(t,to) + i, 

where (^) has been used. This then defines the effective residual interaction 

Mm "■''>. (99) 

t^oc dt ^ ' 

In the lattice simulation the eigenvalues of ^j^ii^to) at asymptotic times 
(within lattice limits) select the low-energy part of the excitation spectrum 
of the two-body system. This will help to suppress intrinsic excitations 
of the composite particles and thus yields an effective residual interaction 
between the mesons (or hadrons) in their respective ground states. How 
well the excitation of intrinsic degrees of freedom is suppressed depends on 
the lattice size and on the design of the meson operators. 
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4.3 Lattice Symmetries 

Lattice symmetries should be utilized to reduce the size of the correla- 
tion matrices C^^^ and C*-^\ For an L'^ x T lattice the action is (usu- 
ally) invariant under the group 0{d,'Z) of discrete transformations of the 
cubic sublattice. The representation theory of these groups is standard 
[ p9| . For example, if c? = 2, there are five irreducible representations 
r = Ai, A2, Bi, B2, E, using a common nomenclature, with dimensionali- 
ties A'^r = 1)1; 1)1; 2 respectively. For d = 3 there are, again in common 
nomenclature, F = Af , , ,T,^ with respective dimensionalities 

A^r = 1)1)2,3,3 . 

Given a fixed, discrete, lattice momentum p' the application of group 
transformations g £ 0{d,'L) generates a set of momenta p = Ogp' that all 
have the same length p = \p\. These transformations define a representation 
of 0{d,'Z) with basis vectors \p > which is in general reducible. Let 

\p>=Y.^\{T,p)e><{T,p)e\p>, (100) 
r e 

where | (F,p)e > denote a set of basis vectors, e = 1 . . ., of the subspace that 
belongs to F. Those can be constructed with group theoretical techniques 
[ p9[ . A simple example, which will suffice for our present purposes, is given 
by choosing p = (p, 0), assuming d = 2 and pL/27r G N, and the trivial 
one-dimensional representation F = ^1, thus 

\{Ai,p)l >= i (|(+p,0) > +1(0, +p) > +|(-p,0) > +1(0, -p) >) . (101) 

Since both the full and the free correlators C^^^ and C^'^\ respectively, 
commute with all group operations g £ 0{d, Z) there exist reduced matrices, 

say C^^'^^ and C^^'^\ within each irreducible representation F such that 
(Schur's lemma) 

< {T,p)e\C<^\t,to)\{T',q)e' >= 5rr'S,,'Cj,f\t,to) (102) 

and similarly for C^'^'^\ Above we identify Cp-g- with < p\C\q >. 

Similar group theoretical considerations are of course also applicable on 
a coordinate space |r > lattice. 

An example from the QED2_|_i model is shown in Fig. The set {W} 
of energy levels, with labels n = 1 . . . 6, belongs to the free system. Those 
were extracted from the asymptotic time behavior of C^^^(t, io)) see (|59|) , 
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and are degenerate (in general) for momenta of the same magnitude but 
different directions. The degeneracy is lifted by the residual interaction, the 
set {W} being extracted from Cp'^q^\t,to). Since the angular momentum 
content of the various representations F = Ai, A2, Bi, B2 is different (Ai 
mostly contains i = 0, Bi and B2 contain i = 2, etc) the hyperfine splitting 
can be attributed to a spin-orbit like force. 

1.7 



1.6 



1.5 



1.4 



11=5 

n-4 



1.3 



1.2 



1.1 



1.0 



W 



0.9 



Figure 9: Hyperfine splitting due to the residual interaction for the QED2_|_i 
model [^0|. The labels for {W} refer to irreducible representations of the 
lattice symmetry group 0(2,Z). 



4.4 Truncated Momentum Basis 

30h that the reduced matrix elements of the free 



It is obvious from (^) and 
correlator have the form 



ci^ \t,to) 



-pq 



.(r) 



it,to) 



(103) 



(V) 

The functions Cp (t,to) linear combinations of 2-point correlator ele- 
ments C*^^/ with momenta p and p' that differ only by their directions. 
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The eigenvalues of C^'^'^\t,tQ) describe the time evolution of the eigen- 
states of a free meson- meson Hamiltonian, say TCq. Since TCo relates to a 
composite system, it is useful to think of its degrees of freedom in terms of 
the relative motion of the two mesons (as clusters) and the intrinsic motion 
of the individual mesons. Assume, for the purpose of discussion, that the 
internal excitation energies are large (hard modes) compared to the relative 
excitation energies (soft modes). In principle, the separation of soft and 
hard modes may be enforced by choosing a big enough L'^ x T lattice: 

• Large L will yield a dense momentum spectrum p oc L^^ . 

• Large T will permit suppression of intrinsic excitations, at large t. 

In practice there are of course numerical limitations, however, one may check 
how well those conditions are satisfied. We assume that d^''^^ and C'^^'^^ 
are available on a truncated set of lattice momenta and both have size 
Naturally this will allow us to study features of the effective residual inter- 
action in the long and intermediate range. 

— (4;r) 

The advantage of working in momentum space is that C ' is already 



diagonal, see (|103| ). Thus the effective correlator (|95D, (|98|) simply becomes 

(Mo) = ,n ^^"^'^^'if . (104) 

cy{t,to) c\\t,tQ) 

We envision numerical diagonalization, say 

Ar(r) 

^fPitM = vP{p)XP{t,to)vP{q) . (105) 

n=l 

The eigenvectors vP are time independent [^], whereas the eigenvalues 
behave exponentially for asymptotic (~) times, 

XP{t,to) ~ aP exp{-wP{t - to)) , (106) 

or, if periodic boundary conditions across the time extent of the lattice are 
imposed, 

xP{t,to) ~ cosh{wP{t - t,)) . (107) 
Now, extraction of the effective residual interaction as defined in (p9|) be- 



(r) 

comes a trivial matter. Provided that the Vn are orthonormal we have 

iV(r) 

< p\nP\q >= ^ vPip) wP vPiq) . (108) 

n=l 



30 



These are the desired matrix elements of the effective residual interaction 
in the truncated momentum basis. They may be used in a variety of ways. 
For example, in the basis |p > of lattice momenta we have 

<p|7^7|^>= J] J]<f|(^,^5)e><^J|7^f)|g>< (r,g)e|9-'> . (109) 
r e 

This would require computation of the reduced matrix elements for all ir- 
reducible representations F of the cubic lattice symmetry group. However, 
depending on the system studied, it may well be that Ai dominates the 
series. 

The coordinate space matrix elements of the effective residual interaction 
are obtained from (|109|) by (discrete lattice) Fourier transformation 



< f\ni\s >= L-^ Yl Yl e'P-^''-'^e-''^<''+^ < p- q\ni\p + q> . (110) 



p q 



It is useful to write the matrix element of TCj as a sum of < —qlTi-il + q > 
and a remainder. A reference system where the relative momenta before and 



after a scattering event are itg, respectively, is known as the Breit frame |61]. 
Then, the sum over p in ( |110[ ) gives rise to 6fg and thus to a local potential 
V(f ) for the Breit-frame contribution, and a genuinely nonlocal potential 
W(r, s ) for the remainder 

<r\ni\s>=6^gV{r)+W{r,s). (Ill) 

The potentials are 

V{f)=Y,e-'^'^-''<-q\ni\+Q> (112) 
q 

W{r,s) = L-'^YYe'P<''-^e-'^<^+^ 
p q 

(< p- q\ni\p + q>-< -q\ni\ + q >) . (113) 

Results from an actual simulation are displayed in Fig. The (lattice) 
Fourier transform, see (110), of (p.09|) is shown. An angular momentum 



projection was employed. Most of the £ = partial wave is contained in 
the Ai sector. The corresponding scattering phase shifts, obtained from a 



Schrodinger equation, are shown in Fig. 11. 
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Figure 10: Example of a meson-meson potential derived within the 
momentum-space approach. Local and nonlocal potentials are shown. The 
model is QET)2+i in 2+1 dimensions |6C]. 



4.5 Adiabatic Approximation 

Probing the lattice with a set of coordinate space operators of the type 

Mt) = Yl = Xd{x, t) Xu{x, t) , (114) 



where we have used (28), gives us an alternative window on the effective 
residual interaction. The Fourier transform of the two-meson operator ( pO| ) 
becomes 

Mt) = e-'P-%{t) = YY1 ^r,^-y Mi) Mt) ■ (115) 

P X y 

It corresponds to two composite mesons with relative distance f and total 
momentum zero. 



Following Sees. 3.3 and we consider the free meson-meson correlator 



77(4), 



(116) 



P Q 
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Figure 11: Scattering phase shifts, for the s-wave, obtained from the effective 
interaction depicted in Fig. The four data points (filled circles) were 
computed in Ref. |^^, using Liischer's method. The same set of phase shifts 
are also shown modulo 180° (open circles). The diamonds on the abscissa 
correspond to the values of = (27r/-L)^/c^ for an L = 24 lattice. 



It may also be written, using (p9D, as 

Cfsit, to) = (e'^^<''~^ + e*P-(^+^^)) \cpit, to)P . (117) 



Here Cp-(t,to) is the correlation function of a single-meson interpolating op- 
erator with momentum p, see (|57|). Its (exponential) time behavior is de- 
termined by the total energy E{p) = \J m? + p"^, assuming continuum dis- 
persion for the moment. In the non-relativistic limit, E{p) — > m + p^ /2m as 
m — > oo, the correlator thus becomes independent of p and may be replaced 



with Cp'=o(i;^o) in ( 117 ). This leads to the approximation 

Cfl{t, to) « L\6r-^g+Sr-^g)\Cp-=o{t, to)|' • (118) 

In coordinate space the free correlator is diagonal in the m — > oo limit, 
whereas in momentum space diagonality is exact. 
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In terms of the fermion propagator we have, see (^) and 
Cp=o(t,io) = j;^|G(xt,yto)|') = (L"' J] |G(ft, xoto)|') • (119) 

X y X 

This is just the time correlation function for a single meson at rest. 

Using lattice symmetry, see Sees. 4.3 and 4.4, the reduced correlator 
matrix has the form 

C^^P{t,to)^Sr,s\c^''Ht,to)\\ (120) 
where the asymptotic time behavior is given by 

c(r)(t,to):^a(^^e-"^(*-*»), (121) 
at least for T = Ai. The full, interacting, correlator matrix is built from 

c(5(t,to) = (4(i) ^s{to)) - {4it)){^,{to)) . (122) 

Again, the separable term vanishes in our example (flavor assignment). Ex- 
pressmg via Wick's theorem, in terms of the fermion propagator and 
the diagrammatic classification proceeds just like in Sees. O and 3^. Al- 
ternatively, we may simply Fourier transform (^), 

d;^it,to) = L-'Y.llT.T.^r,x,-^^5s,y^-y, (123) 

XI X2 yi y2 
(( G*{x2t, mto) G{x2t, mto) G*{xit, yito) G{xit, yito) 
+ G* (fit, y-ito) G{xit, y2to) G* (f 2^, yi^o) G{x2t, yito) 

- G*(x2t,yito) G{x2t,y2to) G*{xit,y2to) G{xit,yito) 

- G{x2t,yito)G*{x2t,y2to)G{xit,y2to)G*{xit,yito))) . 

In keeping with the m — > oo limit it is reasonable to assume that the rela- 
tive distance between the (heavy) mesons will not change much during the 
propagation from to to t. After all, the meson propagation is described by 
( |119| ). This leads us to neglect the off-diagonal elements of C^^\ More 
precisely, since it costs no energy to change the system according to a trans- 
formation that belongs to the lattice symmetry group 0((i, Z), the correct 
approximation is 

Ci^f\t,to)^drsGi^t''^\t,to). (124) 

Now, the eigenvalues of C^'^'^\t,to) are just the diagonal elements. Those 
behave exponentially for asymptotic times. 

C(t^r)(t,to)^4^)e-^'-''(*-*o). (125) 
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The effective correlator ( |104D is also diagonal in this approximation 

4fHt,to)^6rs^if\t,to), (126) 

with 



e:(^;^)(Mo)-^e-K'-H(-o). 

According to our definition ( |99|) the efi'ective residual interaction is 

< r\nP\s >r 



6rsV'^^'> (r) , with V^^^ (r ) = wP - 2m . 



(127) 



(128) 



The above line of arguments is somewhat similar in spirit to the Born- 
Oppenheimer approximation often used with systems of two atoms or mo- 
lecules [^] . The time scale for the dynamics of the interaction is much faster 
than the time scale for the dynamics of the motion for the atoms. Borrowing 
a thermodynamics term, we speak of an adiabatic approximation. In our 
case we simply compute the total energy of the two-meson system at fixed 
(static) relative distance r and interpret the difference to the ground state 



energy as a potential ||64|. 

An example from numerical work on the QED2_(_i model [65| is shown in 
Fig. |l^. The left part displays the static potential V*-'"^ according to ( |128D . 
The distance r = is forbidden by the Pauli exclusion principle. This is the 
reason for the increase (repulsive core) of the potential as r — > 0. 

The right part of Fig. [l^ shows a potential computed slightly differently, 
for the sake of comparison. The perturbative version of (|T2^) leads to 



V(r)(r 



t-to 



r(4;r) 



1 



{t,to) 



aP\-a(^)\ 



(129) 



In this way the r-dependence of a)- enters the analysis. The results in 
Fig. ^ are similar. They shed some light on the validity of the perturbative 
definition of Tij put forward in Sec. 4.1. 

Of course, only the local part of the effective residual interaction is ac- 
cessible in this way. Note, however, that a local potential 



V(f ) = ^ ^ < ^|(r, r)e > 



(130) 



compared to (109), may not be the same as the approximate local potential 



of (112) even if the untruncated correlation matrices are used. The rea- 
son is that diagonal elements of the non-local part of the effective residual 
interaction may contribute to the adiabatic local potential. 
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Figure 12: Static meson-meson potential for the model according 

to ( |128| ), on the left. On the right a somewhat different approximation 
according to (129) serves to check the perturbative definition of TLj. 



The adiabatic approximation hinges on the assumption of heavy (m — > 
cxd) partners. This should be reasonable if at least one of the quarks in each 
hadron is heavy, say an s-quark. However, in the chiral limit (m^r 0), 
using Wilson fermions, the adiabatic approximation is expected to fail if 
pseudoscalar mesons are involved. In this case momentum space methods 



like in Sec. 4.4 are preferable. Finally, the quality of the adiabatic ap- 
proximation can of course be tested numerically by computing off-diagonal 
elements of the correlation matrices. 

4.6 Analysis on a Periodic Lattice 

Periodic boundary conditions across the space extent of the lattice are a 
common choice. The potentials V(r) and >V(r, s), see (|112|) and ( |113D , 
extracted from the lattice simulation, in one or the other way, will reflect 
those conditions. In particular, the maximal usable relative distance r is 
L/2. 

Suppose it is desired to make a fit to the lattice local potential with a 
class of functions V^'^\r) depending on a set of parameters a. An example 
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IS 

(f) = ai ^ + CO . (131) 

This class (see Fig. |l^) has enough flexibihty to match heuristic features of 
the hadronic interaction, hke short-range repulsion or attraction as the case 
may be, and a long range Yukawa form 

y(")(f) ^ -^i^^^ + ao as r^oo. (132) 
as r 

The periodic extension, defined as 

y^°) (f ) = ^ y (f - nL) , (133) 

n 

where n are vectors with d whole number components, can then be used to 
fit the lattice potentials with 

^'(") = ^ E - yt\r)\'^ . (134) 

f 

Here the sum runs over the set of A'^ lattice distances f for which V(f ) has 
been computed. 

The idea behind using a periodic extension, like ( 133 ), is to account 
for the fact that the lattice hadrons are interacting with their replicas in 
adjacent copies of the periodic lattice. With minimized, we would take 
\/{")(f) and W^f^\f, s) obtained from a similar fit, as the results of the 
simulation. 



5 Current State in QCD in 3+1 Dimensions 

At the time of this writing, lattice work on the subject of hadronic inter- 
action is in an exploratory phase. Nevertheless, it is useful to convey what 
has been done through some selected examples. We will not include in this 
section results other than QCD3_,_]^, for example 0(4) symmetric field the- 
ories and four-fermion models |6^, |6^, models exploring static quark 
geometry^ [7^, studies with an SU(2) gauge group ||7l[|, and work using 
the hamiltonian lattice formulation j72[ |. 

^The contribution by A. M. Green is devoted to this subject. 
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Figure 13: Example of V^"^ as defined in (131). The parameters a are 
arbitrarily chosen, for the purpuse of illustration. 



5.1 Scattering Lengths for tt and N Systems 

To date the most elaborate application within QCD3^^ has been via Lii- 
scher's formula (|9|) by the CP-PACS collaboration aiming at scattering 
lengths for the vr-vr, vr-N, K-N, K-N, and N-N systems [||, [l|, 



24, 25, 26 1 . The collaboration has applied considerable resources to 
this problem, recently using lattices as large as x T = 32^ x 60 |2|]. 



However, for most of their earlier explorations |21] lattices typically of size 
12^^ X 20 with lattice lattice constants in the range of a ~ 0.15 - 0.30 fm 
were used. 

Extracting the s-wave scattering length oq from Liischer's formula 



Whih2 - (mhi + ruh^) = - 



-3 [1 + Cl— + C2(— ) +o{L 



(135) 

where hi and /i2 denote hadrons 1 and 2, see (^), requires a standard mass 
calculation on the lattice. One therefore needs to compute correlation func- 
tions of suitable operators, say (phiit) and (j)h2{i)- For illustration, some 
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examples using Wilson fermions are 

'Pw+it) = -'^ipdAix,t)-f5ipuA{x,t) (136) 

X 

(t>7To{t) = --^^{ipuA{x,t)j5iJuA{x,t) -'lpdA{x,t)j5'll;dA{x,t)) (137) 

X 

4>K+(.t) = '^i^sA{x,t)j5^uA{x,t) , (138) 

X 

etc, where A = 1,2,3 represents SU(3) color. Only zero-momentum in- 
terpolating fields are needed here. Kogut-Susskind, or staggered, fermions 
p7| , p8| , |49| ] were also used by the CP-PACS collaboration. The construc- 
tion of operators with specific hadron quantum numbers is technically more 
complicated, however [[7^]. The authors of Ref. [21| chose to construct cor- 
relation functions where sources and sinks between hi and /12 are one time 
slice apart. For the one-hadron systems this means 

C,,(t,0) = (0/.,(t)<(O)) -(</.,,(*)) (4^(0)) (139) 
Ch,{t + 1,1) = (0,, (t + 1)4^(1)) (t + 1)) (4^(1)) (140) 



where ~ indicates the asymptotic behavior. The two-hadron interpolating 
fields are constructed from linear combinations of products of one-hadron 
operators, respecting the one-time-slice offset. Denoting those by ^hih2 (^i ^+ 
1) illustrational examples in the tt-tt system with isospin / are 

^ii='''^='\t,t + l) = l=(^^^{t)cp^_(t + l) (141) 

-0,0 {t)(t>^o (t + 1) + 0,- (*)(/>,+ (t + 1)) 
^ii=^''^=^\t,t + l) = <A,+ (*)</.,+ (t+1). (142) 

In the s-wave I = 1 is not allowed (Bose symmetry). The corresponding 
two-hadron system correlator is 

Ch,hAt + 'i-,t,l,0) = (143) 
{^h,h,{t,t + !)</,, (0, 1)) - {^h,h,{t,t + !))(</,, (0, 1)) . 

The measurement of the s-wave scattering length oq via Liischer's formula 
(|135| ) requires the energy difference Wh^h2 ~ ("^-fti + "i/ij), which can be ex- 
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tracted from the asymptotic time behavior of the ratio of correlation func- 
tions 

C,,(t,0)C,,(t + l,l) - a,,a,,' ■ ^ > 

Baryon interpolating field operators are more complicated. For the pro- 
ton (A^^) and neutron (N^) standard choices for zero-momentum operators 
are 



')N+ 



9^0 



(t) = ^ eABC (jPuA (x, t)C^dB {x, t)^ IpuB {x, t) 

X 

- (^'tjjJ^{x,t)CipuBix,tj^ ipuB{x,t) (145) 

(i) = ^(-ABC (i'uAix,t)C^dBix,t)^ 1pdBix,t) 

i^dAix,t)C^uBix,t)) iJdBix,t) , (146) 



where C means charge conjugation and A, B, C are color indices. Operators 
for the "^Si and ^Sq channels of the N-N system are constructed in Ref. [^] 

«>35,(t,t+l) = -^(</.jv+(t)0^o(t + l)-</.^o(t)(/>;V+(t+l)) (147) 
cl>i5„(t,i + l) = ct)N+{t)(t>N+{t+l). (148) 



Again, ratios like ( |144| ) allow the extraction of the desired mass shifts. 

Computation of the above time correlation functions requires matrix 
elements of the quark propagator G = D~^, where D{x,y) denotes the 
Dirac, or fermion, matrix and x = (x, t), are space-time lattice sites. Since 
a complete solution of 

^ D{x', x") G{x",x) = 1 6,,^^ , (149) 

x" 

where 1 means the unit matrix in color-Dirac space, is not feasible in practice 
source techniques are typically used to obtain at least some columns of G. 
The CP-PACS collaboration has employed so called wall sources. These are 
defined by placing a point source of value 1 on each spatial site of a certain 
time slice t, specifically 

D{x't', x"t") Gt{x"t") = 1 hs't'Wt) = 1 ^t't . (150) 

x"t" X 
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Note that the uniform nature of the wall source implies that the solution 
vector Gt must be independent of x' . Indeed, employing ( 149 ) it is obvious 
that 

Gt{x"t") = Y G{x"t", xt) . (151) 



Loosely speaking, since the spatial source is structureless quark propaga- 
tion to x" proceeds in 'the same way' regardless of the spatial point of 
origin. However, this is not a gauge invariant mechanism. Under a gauge 
transformation D{x,y) = U^x)D{x,y)U{y), etc, the wall source becomes 
lU{x't') 6ti Thus, in contrast to G{x"t" ,xt), the (numerical) solution 
Gt{x"t") of ( |150| ) does not transform covariantly. Gauge dependent noise 
in hadronic correlation functions caused by the wall sources is specifically 
created at time slices where two wall sources or one wall source and one 



sink are placed. This can be seen from Fig. |lj which shows a diagramatic 
analysis of the correlators for the vr-vr, vr-N, and N-N systems in terms 
of quark propagators and the corresponding wall sources and sinks. Wall 



t, - 
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(ii) 
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Figure 14: Diagrams contributing to hadron-hadron 4-point functions, from 
Ref. |21|. Arrowed lines represent quark or antiquark propagators. Wall 
sources are depicted by short bars, and circles represent sinks for hadron 
operators. The diagrams in (i), (ii), and (iii) are for the vr-vr, vr-N, and N-N 
systems, respectively. The decompositions (|153|) and ( 154 ), for example, 
refer to the labels D,C,R,V of diagrams (i). 
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sources are depicted as short lines and sinks as circles. The time slices refer 
to ti = 0,t2 = 1, = t, i4 = t + 1. For example the vr-vr rectangular diagram 



in Fig. 14 (i) corresponds to 



C^(t4,t3,i2,tl) = Yl (ReTV[Gl^(x2t2)Gt,(x2t2)Gljf 3*3)^^1 (f 3^3)]) ■ 



(152) 

The paths involved in the computation of are not all closed, they are 
interupted at time slices where a (nonlocal) wall source meets a (local) sink, 
or two wall sources meet. This happens, for example, at ti and t^, see (|152|) 



and Fig. |14. 

In order to deal with the gauge noise problem two strategies may be 
employed: The first one is gauge fixing, the other one is to rely on gauge 
fluctuations to cancel out in the final correlation function. For the latter the 
gauge noise is expected to decrease as ~ L^/^ for sufficiently large L |21]. 



Correlator ratios ( 144 ) for various channels can be expressed in terms of 
linear combinations of certain diagrams which are referred to by labels D, 



C, R, V, etc, in Fig. 14. For example, in the Wilson fermion scheme 



RDit) + lRO{t)-m''{t) + lE^{t) 
R^{t)-R^{t). 



(153) 
(154) 



For staggered fermions the corresponding decomposition is slightly different. 

Typical results for a ratio function (144) are shown in Fig. |l^ for the 
7 = and 1 = 2 channels of the vr-vr system. A selection of scattering 



lengths, again for the tt-tt system, is displayed in Fig. |16|. In Fig. |17| a com- 
parison of scattering lengths for the systems tt-tt, tt-N, and N-N is made. 
Table |l| contains the main results of the lattice simulation and current al- 



gebra predictions. A chiral extrapolation has not been made in |21]. Thus, 
perhaps not surprisingly, the lattice results agree best with the current al- 
gebra predictions that use the hadron masses supplied by the actual lattice 
simulation. The N-N system is special because of its large size. Realistic lat- 
tice results for the scattering length cannot be expected with contemporary 
resources. 
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Figure 15: Ratios (153) and ( |154| ) for the vr-vr four-point function, from 
Ref. ||2l|. This example is for Wilson fermions at (3 = 5.7, k = 0.164, and 
wall sources without gauge fixing. The lattice size is 12'^ x 20. 




0.0 0.2 0.4 0.6 



Figure 16: Comparison of various scattering length results in the / = 
and 1 = 2 channels for the tt-tt system, with units as given in the picture. 
Solid and open plot symbols refer to Wilson and Kogut-Sussking fermions, 
respectively. Results are from Ref. [21|, except for the squares and diamonds, 
which are from Refs. ||l^, Dotted lines are current algebra predictions. 
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Figure 17: Comparison of scattering lengths for the vr-vr, vr-N, and N-N 
systems in physical units oq [fm] versus the n-p mass ratio. These are 
(quenched) results for Wilson fermions from Ref. [21 1 at lattice constant 
a = 0.137(2) fm. 



Table 1: Selected results from Ref. |21|. Scattering lengths ao for the s-wave 
in various channels are shown in physical units [fm], for the lattice (lat), 
experiment (exp), and current algebra (cua). The lattice results are for 
Wilson fermions on a 12'^ x 20 lattice with P = 5.7. Light-quark masses are 
large as indicated by m.,^/mp = 0.74 and rriN /nip = 1.57. Experimental and 
current algebra results are compiled from Refs. [74, 75| and [^], respectively. 
In the last column the masses from the lattice simulation enter the current 
algebra predictions. 











LAT 


EXP 


CUA 


cua(lat) 


TT- 


-TT 


/ = 





+3.02(17) 


+0.37(7) 


+0.222 


+3.47(5) 






/ = 


2 


-0.924(40) 


-0.040(17) 


-0.0635 


-0.993(16) 


TT- 


-N 


/ = 


1/2 


+3.04(66) 


+0.245(4) 


+0.221 


+2.701(41) 






/ = 


3/2 


-1.10(20) 


-0.143(6) 


-0.111 


-1.350(20) 


N- 


-N 






+10.8(9) 


-5.432(5) 










So 




+8.0(1.1) 


+20.1(4) 






K- 


-N 


I = 





+0.55(47) 


-0.0075 












I = 


1 


-1.56(13) 


-0.225 


-0.399 


-2.701(41) 


K- 


-N 


I = 





+4.64(37) 


-1.16 + i0.49 


+0.598 


+4.051(61) 






I = 


1 


+2.63(64) 


+0.17 + f 0.41 


+0.199 


+1.350(20) 
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5.2 Static N -N and N - N Potentials 



Beyond the zero-momentum limit (scattering lengths) systems of hadrons 
with at least one heavy quark each pose the least technical difficulties. Going 
to the extreme, an infinitely heavy quark may be interpreted as a static color 
source. This idea goes back to Wilson's original work [^] and played a role 
in demonstrating confinement. 

The time evolution operator of a static color source (a quark) located at 
fixed X is given by the product of link variables U4{xt) € SU{3) along a line 
U ^ ti + 1 ^ . . . — > in the time direction, ^ = 4. For a lattice x T 
with periodic boundary conditions in the time direction the line is closed, 
ti = tf. The trace of such an operator 

1 ^ 

Q{x) = -TvHu^m (155) 

t=i 

is gauge invariant and is known as a Polyakov loop. It represents a static 
quark at x, and Q^{x) represents a static antiquark at x. In the finite- 
temperature formalism^ the expectation value |]7^ , 79, 64] 

^(Q(f + x)Q\x)) = e-^^QQ^''^ (156) 

X 

is related to the free energy FQQ{f), loosely (but imprecisely) referred to 
as the heavy quark-antiquark potential. Adopting this interpretation it is 
interesting to compute the free energies of two clusters of two or three quarks 



q, 81, Consider 



Oix) = V p{yiy2m)Q{x + yi)Q{x + y2)Q{x + m) , (157) 



where p is some function symmetric in yiybys, for example a Gaussian or 
the characteristic function of a sphere with radius R, which serves as an ad- 
hoc model for the three-quark cluster. The free energy of two three-quark 
clusters at separation r then is 

F{r) = -^lnL-^'^{0{f+x)0{x)). (158) 



An example of F{r) can be seen in Fig. 18, There is attraction in the 
overlap region of the clusters, somewhat reflecting the shape of p. The 
cluster interaction possesses no dynamics in this construction. 

^The temperature is the reciprocal of the time extent of the lattice, 1/Ta. The path 
integral over fields becomes periodic in euclidean time |77[. 
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Figure 18: Interaction energies of hadron-hadron systems versus the relative 
distance ||8^, ^ in the finite-temperature formahsm. Hadronic structure is 
represented by a static cluster. The radius of the cluster is i? = 2 in lattice 
units. 

5.3 Heavy-Light Meson-Meson Systems ^ 

A step towards a more realistic calculation is to make some of the quarks 
dynamical. In a system of two mesons which consist of a heavy and a light 
quark each, we may still approximate the heavy quark by a Polyakov loop, 
as above. The role of the static quarks then is to localize the mesons. This 



situation suggests using the coordinate-space approach of Sec. 4.5 
Thus we take, say in the staggered scheme, 

(pxit) = Xh{x, t) Xu{x, t) , (159) 

where h is the heavy flavor i(/i = s,c, ...). The two-mesorJ interpolating 
field is defined just like ( |115D . However, the heavy- flavor propagator. 
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is taken in the limit rrih oo. This is achieved by using the hopping 

parameter expansion [S3, 77 1 and keeping the leading term only. We thus 
use 

GW(xt,yto) = ri ^^(i'r) . (160) 

r=l 

The product of link variables becomes proportional to the Polyakov loop 
( [155| ) for t-to = T. With G'-^^ diagonal in the space sites the 2-point and 
4-point correlators become 

C(2)(t,to) = {L-'^Y.^BAi^t^^to)GBAixt,xto)) (161) 

X 

cf\t,tQ) = {L-'^Yl[G^l;]^{xt,xto)GBA{xt,xto) (162) 

X 

X G^^*(x — ft, X — fto) Gdc{x — ft, X — f to) 



G^ji {xt, xto) Gbc{x - ft, xto) 

^DC 



X G^^*(f - ft, X - fto) GDA{xt,x - fto)]) 



Sums over color indices A,B,C,D are understood. The 4-point correlator 
( |162| ) is a sum of two terms which are illustrated in Fig. |l^. The light 
quarks can be exchanged between the mesons. Thus, besides gluon effects, 
also quark dynamics is included in the effective residual interaction. 




Figure 19: Diagrams for the heavy-light meson-meson system. The thick 
and thin lines represent propagators for the heavy (static) and the light 
(dynamical) quark, respectively. 

For distances f / parallel to the coordinate axes, the Ai sector pro- 
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jection, see Sec. 4^, is similar to (101), thus 

r,\r\=r 

According to ( |127| ) the effective correlator 

Cr^'^^\t, to) 



(163) 



V(-*i)(r)(t-to) 



(164) 



|C(2)(t,to)|2 

serves to extract the (adiabatic) potential V^'^'^\r). Examples from an ac- 



tual simulation |33] are shown in Fig. pQ. 
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Figure 20: Meson-meson energies W{r) = y(^i)(r), see ( |164| ), for several 
light-quark masses of staggered fermions. Energies at distance r/a = 3 are 
degenerate, |f| = |(2,2, 1)| = (3, 0,0)|. Those data points are slightly shifted 
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In the simple case of a one-hadron local operator it is possible to express 
the corresponding time correlation function in terms of propagator matrix 
elements from a source located at only one fixed lattice site, say x = (xq^o)- 
This is done by using translational invariance of the lattice action in the man- 
ner outlined in Sec. 3.2 leading to (40), for example. The cost of computing 
this is greatly reduced because only one column G{xt, xoto) of the inverse 
fermion matrix (per color-Dirac index) is needed. In the case of a two-hadron 
system translational invariance might still be utilized to replace some prop- 
agator elements with G{xt, xoto), for example via xi xi + yi — xq in (|4^ . 
In the latter case elements G{xt, yto) where x and y run separately over all 
spatial lattice sites still remain. Those are known as all-to-all propagators. 
Their computation poses an almost unavoidable problem in multi-hadron 
systems. The wall sources defined through ( |150D seemingly circumvent this 
difficulty, however, only at the cost of violating gauge invariance of the cor- 
relation functions, see Sec. 5^. In a heavy-light hadron-hadron system, like 
( |162| ) for example, the needed minimum number of sources is the number of 
relative spatial distances distances f that are probed (including those related 
by 0(3, Z) symmetry). 

The use of stochastic methods for dealing with the all-to-all propagator 
problem was already discussed in some detail. Even in cases where transla- 
tional invariance can be utilized it may be advantageous to employ random 
sources instead. The reason is that the additional (large) sums over lattice 
sites work to reduce statistical errors in the correlation functions. Besides, 
it may also be the only computationally viable strategy. In a typical appli- 
cation a random source is chosen to be nonzero on one time slice only, say 
to, as in ( |65[ ) and (|6^). In this case propagation from any lattice site to any 
other lattice site within the same time slice is computed using the stochastic 
estimator, whereas propagation between lattice sites with different times, 
such as t and Iq, is computed with machine precision. An important con- 
sequence is that the errors for propagation in space and time directions are 
different. While the latter is as small as the matrix inversion algorithm per- 
mits (e.g. a conjugate gradient) the former is limited by the variance of the 
random source probability distribution and, in particular, is independent of 
the distance between the lattice sites. This is a serious drawback, it means 
that stochastic estimation is futile in computing the usual (hadronic) time 
correlation functions since their exponentially decreasing signal would be 
concealed by the essentially constant error of the stochastic estimator. 

This problem can be dealt with to some extent by a method called maxi- 
mum variance reduction [35|, see also |86, 87 1. The idea is to enclose disjoint 
regions of lattice sites by fixed-source boundaries in order to increase sample 
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sizes. One such region may consist of a single point ||8^ . Another, somewhat 
related, example is ii = {{x, tQ)\x £ L^}, i.e. the spatial volume at time slice 
to- An estimator for propagation between any two points in disjoint regions 
is then variance reduced. Specifically, starting from 

Z = J [##*] exp{-^^^A^) , (165) 

where (j) is an auxiliary scalar complex lattice field and A is a hermitean 
positive definite matrix, one may generate stochastically independent en- 
sembles of fields {(f)} by way of common Monte Carlo methods with the 
gaussian probability distribution implied in ( |165| ). Writing i = Cfix, etc, for 
the color-Dirac-site indices it is then straightforward to show that 

A^' = I / [##1 eM-\(t>^M) ^j4>i = (0*0i> , (166) 

where in this context (...) means the average over the ensemble {(j)}- The 
above equations become useful for the choice 

A = Q^Q, (167) 

with Q = Q [U] being the fermion matrix for a given gauge field configuration 
U. In the Wilson scheme, for example, it has the form Q[U] = 1 — kD[U]. 
Quark propagator elements Gij can then be estimated from the ensemble 
averages 

G^/^' Qr^'^{{Qcl,)*<P,). (168) 

If the variance of the above estimator is aij the statistical error is Uijj -^N^ 
for an ensemble of size A*"^.]^ However, the crucial point here is that the 
variance cjjj, for i = Cfix and j = Buy, is essentially independent of the 
space-time distance \x — y\, say a^j ~ a while a is of the order of one in- 
herited from its gaussian ancestry. This effect spoils using the stochastic 
estimator for the purpose of computing exponentially decreasing time cor- 
relation functions. 

To alleviate the problem the authors of Ref. |8^] proceed to consider a 
subset, or region, of contiguous lattice sites. Describe this region by sets of 
sites R and S, with i? n = 0, where R contains all entirely interior points 

^In fact, since a new estimator is computed for each of the gauge field configurations in 
the lattice simulation, the resulting error of a hadronic time correlation function is much 
smaller. For this reason a value of A^^ that is about 1/10 - 1/20 of the number of gauge 
field configurations is usually sufficient. 
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(having 8 neighbors £ RU S) and S is the boundary (having 1 through 7 
neighbors G i?U 5"). In the process of generating ensembles {i;^} from ( |165D , 
consider a certain field (p. Then, the field conponents on the boundary 5 
are kept fixed 



si =^ (pi for / G 5 , (169) 



and instead of (|165|) one uses 



Z = JjdcPdcP*] exp (^-^{4Aij4>, + 4Ausi + slAij^,)^ , (170) 

where i,j € R and I £ S is understood for the index sum ranges and the 
notations A and A are introduced for the corresponding blocks of A. The 
integral runs over the field components located on R. Now the i?-averaged 
random source field 

Vk = ^ jjd4> d4>*] exp (^-^{^lAij(l)j + (t)\AuSi + s\Uj(pj)^ 4>k , (171) 

with k £ R, is called the variance reduced estimator for (p. Evaluation of 
the gaussian integral in ( |171| ) gives 

Vk = -Ag'AiiSi. (172) 

Recall that in ( |168| ) single-site components of the random source field 
(p are employed to estimate the quark propagator. Replacing those with 
site-averaged variance reduced estimators v essentially amounts to using 
extended sources which of course reduce fluctuations. This is the main idea. 
To obtain a variance reduced estimator for quark propagator elements two 
disjoint regions R and R' are needed. Then, within the above framework, 



one may use (168) to show that 

Gk'k - {{Qv)l v[,) for k£R, k' eR' , (173) 
where v and v' are the variance reduced sources corresponding to the regions 



R and R' respectively |87]. Clearly, only propagator elements between any 
site in R and any site in R' can obtained in this way. Thus the all-to- 
all propagator problem is only partially solved. In order to maximize the 
variance reduction effect, R and R' should be equal in size and together with 
S and S' cover the entire lattice. This suggests the choices of all lattice sites 
with < i < T/2 and T/2 < t < T for R and R' , respectively. Further 
computational and technical details for the practitioner are given in Ref. |85]. 



51 



Michael and Pennanen have utihzed maximal variance reduced estima- 



tors in their work on heavy- light meson- meson systems |34] using improved 
Wilson fermions with a Sheikholeslami-Wohlert action [^] . Their main mo- 
tivation was to explore the possibility of binding heavy-light systems. Since 
the kinetic energy of the quarks has a repulsive effect at shorter relative 
distances two hadrons containing a heavy quark each should have a better 
chance of forming a bound state. Concerning experiment, the simulation 
is expected to describe mesons with one (heavy) b-quark flavor. Using the 
static limit for the heavy flavor has the consequence that heavy-light pseu- 
doscalar and vector mesons (B and B* respectively) are degenerate. How- 
ever, experimentally their mass difference is less that 0.9% of the ground- 
state B-meson mass |89], so the static approximation is quite reasonable. 



The notation B is used in in Ref. |34] to refer to the degenerate heavy- light 



mesons. For the light quarks u and d flavors are employed, the correspond- 
ing Wilson hopping parameter k puts their masses in the neighborhood of 
the strange quark mass. 

Thus a B-B system with relative distance r between the (static) quarks 
is classified in terms of the total spin Sq and isospin Iq of the light quarks, 
while imposing the usual particle exchange symmetries through appropriate 
heavy quark total spin Sb and spatial symmetries. The latter states are 
degenerate in the static quark limit. In the limit of an isotropic spatial 
wavefunction (L=0) there are only four degenerate ground state levels of 
the B-B system characterized by /g = 0, 1 and Sq = 0, 1. The situation is 
illustrated in Table |2|. 



Table 2: Allowed B-B states with L = according to |34|. Total isospin 
and spin of the light-quark subsystem are denoted by Iq and Sq, respectively. 
The total spin Sb of the heavy quarks only distinguishes degenerate states. 
Also shown are the physical heavy-light systems and their coupling to 
these states. 
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Selected results of the simulation |34] are shown in Fig. 21 exhibiting 
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the energies of the B~B system at fixed separation r of the static 6-quarks. 
The authors suggest that there is evidence for deep binding at small r with 
a light di-quark configuration having {Iq, Sg) = (0,0) and (1, 1) respectively. 
This binding energy is 400 - 200 MeV at r = but is very short-ranged. It 
is essentially a gluonic effect and is rather insensitive to the light quark mass 



as shown by studies of static baryons with varying light quark masses |85]. 
Possibly, a configuration where the color of two heavy quarks combines to a 
color state symmetric under particle exchange (sextet) provides the interac- 
tion mechanism here. This situation is particular to relative distance r = 
since two-quark color states other than local singlets are not suppressed by 
confinement dynamics. 
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Figure 21: Results for the binding energy E[B-B) — 2E{B) between two B 
mesons with light quarks in (Ig, Sq) = (1, 1), (0, 0), (1, 0), (0, 1) at separation 
R (called r throughout this chapter) in units of Rq ~ 0.5 fm. The light 
quark mass used corresponds to strange quarks. Results at different spatial 
lattice sizes are displaced in R for legibility |^ . 
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At larger r, around 0.5 fm, evidence is seen for weak binding when 
the hght quarks are in the {Ig,Sg) = (0,1) and (1,0) states. This can be 
related to meson exchange and one finds evidence of an interaction in the 
spin-dependent quark-exchange (cross) diagram which is compatible with 
the theoretical contribution from pion exchange in the study. Using lighter, 
more physical, light quark masses below the strange quark mass and, most 
desirably a chiral extrapolation, will be necessary to confirm this interaction 
mechanism. 



5.4 Momentum-Space Work on the vr— tt System 



The calculation of euclidean correlation functions between different hadrons 
in the incoming and outgoing channel, like K— > vrvr decays, has to be treated 
with caution ||9^. The program laid out in Sees. 3^-3^ and has been 
put to the test in QCD3^^ for a tt-tt system of Wilson fermions . The 
formalism is essentially unchanged, except for the appearance of color indices 
on the link variables U^{x) G SU(3) and, replacing staggered with Wilson 
fermions, the appearance of both a color A,B... = 1,2,3 and a Dirac 
H,^ ... = 1,2,3,4 index on the Grassmann quark fields ipfA^ix). In addition, 
there is a rather technical, but important point: 

A technique known as 'smearing' uses spatially extended field operators. 
For example, omitting flavor and Dirac indices for the moment, define the 
quark field recursively 



A'Hxt) 



x,y 



(174) 



with k gN, and the matrix 



KAB{x,y\ 
HAB{x,y] 



^AB 5x,y + OiHAB{x,y) 



(175) 



E \yrnAB{xt)5s,y~rh + Ul^^j^Q{yt)5s^y+ih (176) 



m=l 



connecting (hopping) to next-neighbor sites. With this particular choice of 



the matrix K the procedure is known as Gaussian smearing |91]. The real 
number a and the maximal value for k are considered parameters. The ad- 
vantage of smearing lies in the fact that interpolating operators constructed 
with smeared fields often have a larger overlap with the ground state of the 



hadron in question, see Fig. 22. The signal from the corresponding correla- 



tor is less contaminated by fluctuations from excited states and thus masses 
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can be extracted from shorter euclidean time extents. This leads to better 
statistics and is often crucial for the numerical work. A similar technique, 
known as 'fuzzing' applies to the link variables [^3| but will not be detailed 
here. 



k=0 k=l k=2 



Figure 22: Illustration of the 'smearing' technique for fermion operators. 
Various iterations {k = local, etc) are shown. The thickness of the lines 
indicates the 'multiplicity of use' of a particular link. 

Thus, the equivalent of the one- meson field (|2^) is replaced by 

^t) = ^ {x, t)75V'S t) (177) 

X 

for a TT"*" pseudoscalar meson. As usual, the quark propagator is obtained 
from contractions between the (unsmeared) quark fields 

n 1 

• • • V'/Am (x) ^gBv (y) • • • = • • • GfA^,gBu{x, y) . . . (178) 

■n n 

■■■^fA^,{x)^gBu{y)--- = ■■■-l5,,i't,G}A^,^gBu'{x,y)-f5^^,^> ... . (179) 

The latter equation is a consequence of the identity 

G(y,x) =75Gt(x,y)75. (180) 

The t refers to flavor-color-Dirac indices. The identity is valid if the lattice 
action is CPT invariant ||9^ . Further, if the mass matrix in the Wilson 
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fermion action is flavor diagonal, then also 



GfA^,,gB^. {x,y) = SfgC^ll s^ {x,y) . (181) 

Typically, one would choose equal Wilson hopping parameters k (those de- 
termine the quark masses) for u and d quarks and different ones for the 
strange quark. 

We now give a list of essential things to do, according to the plan of 
Sees. I and |. 

Consider the linear equation 

yy4 Bu 

(182) 

The meaning of the indices are A^B... = 1,2,3 color, /i,z^... = 1,2,3,4 
Dirac, x,y . . . space {d = 3), X4, time, and r = 1 . . . Nji labels the random 
sources R for each source point. A prime ' denotes a source point. There is 
some freedom in choosing the latter. In ( |182D the sources are nonzero on one 
time slice only. The same set of sources is used for different flavors (Wilson 
hopping parameters k), but a separate source is chosen for each color, Dirac 
and time index. The features of the random sources are somewhat linked 



to the choice of the solution algorithm |9S]. Also, Gaussian, Z2, or other 
random sources may be employed. In any case, the random-source average, 
which we approximate numerically as 

(r) « r=l 



must satisfy 



(r> 



A conjugate-gradient algorithm |94|, or a variant, is suitable for solving the 



linear Eq. ( |182| ). An estimator for the propagator matrix elements then is 

<,A,(y"*2/4,xx4) =Y.xl,f^^^-^\yy,)R(^>^^^^>{x) . (185) 
{r} 

For simple hadron-hadron systems it may be sufficient to place sources on 
one time slice X4 only. This is the situation for the TT~^~7r~^ system (/ = 
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2 channel). Matching the notation in the examples discussed in previous 
sections we set X4 = tQ. Putting things together, the correlation matrices 
are then computed through the following sequence of steps: 



rP^^''''''\x) = 5cbR^^''''''\x) (no sum overs) (186) 
i?«(^^-')(x) = Y^Y.^cAix,y)Rt'^^''''^''\y) (187) 



X^ff'^-^'\,t)=xiJ;^'-'\st) (188) 
y A 

X C 
^^W{AMr'*')*(^)^W(B.rt')(^^ (190) 

X C 

Summation over doubly-occuring indices is understood, except indicated 
otherwise. The expressions for the 4-point correlators C^^^ and C^^^ are 
more complicated, but similar. 

Figure ^ shows energy-level diagrams of the non-interacting, from C^^^ , 
and the interacting, from C^'^\ tt-tt systems for two values of the hopping 
parameter, = 5.888 and = 5.972. These are from a set of six 
values which have been used to perform the chiral extrapolation 
(equivalently k — > Kc). The other lattice parameters were L = 9, 
T = 13, and f3 = 6.2 with the next-nearest-neighbor o(a^) tadpole-improved 
action [ pT| , The lattice constant, matched to the string tension, is a = 
0.4 fm or = 500 MeV. An interesting feature of the spectra is that 
high-momentum levels are lowered by the interaction indicating attraction 
at short distances, and some of the lower levels move up indicating some 
repulsion in the system. 



The corresponding approximate local potentials, according to ( 112 ) of 
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Figure 23: Meson-meson energy level splitting due to the residual interac- 
tion. Two different values of the hopping parameter (quark masses) are 
shown [S2|. 



Sec. 4.4, were projected onto the partial wave £ = 

Mr) = d^r-Vir) = M2qr) ^^k) , (193) 



with q 



0...k„ 



Results of the chiral extrapolation plot |p5[ are 



shown in Fig. 24. The Fourier nature of the obtained potentials is apparent. 
A parametric fit, using (|l3Tl),(p^) and ( |l3l , is also shown. Of course short 
distances r are beyond the scope of the truncated momentum basis and 
should be interpreted with caution. Thus only small relative momenta make 
sense when calculating scattering phase shifts with this potential. Results 
using standard (non-relativistic) scattering theory p^ ] with parameterized 
potentials ( 13lD are displayed in Fig. 
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T / a r/a 

Figure 24: Approximate s-wave projected local meson-meson potentials ac- 
cording to (|112| ) and (|193|) . Shown is the result of a chiral extrapolation 
(m^ — > 0) of the potential based on various hopping parameter values; curves 
resulting from the use of 3 to 6 k- values are compared. A parametric fit with 
y(")(r) based on (|l3ll ), using 3 values of is presented pSj. 




0.2 0.4 0.6 0.8 1-0 1.2 1.4 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

M^^ [GeV] M^^ [GcVl 

Figure 25: Scattering phase shifts 5jzQ calculated from the tt-tt potential 
of the lattice QCD simulation (thick lines) |^5|. Influence of momentum 
cutoff /cmax = 3 (left) and /cmax = 4 (right) is shown. Results using the 
classical and relativistic dispersion relation are distinguished by m-,^ and 
m-,r{p), respectively. Errors are represented by the dotted regions. Their 
boundaries (dashed lines) correspond to the phase shifts calculated with 
extremal (bootstrap) potentials. The experimental data are from Ref. |Q. 



59 



5.5 Coordinate-Space Work on the tt— tt System 

We turn to another exploration in coordinate space with staggered fermions 



[ p2[ . The final interaction potential is shown in Figs. 26(a) and (c) for nif = 
100 MeV and ruf = 500 MeV, respectively. The fits of the periodic potential 
were stable for the data with ruf = 100 MeV. The data for mf = 500 MeV 
suffer from very small energy differences which makes a fit unreliable. The 
potentials exhibit attraction at intermediate range and may have a repulsive 
behavior at short distances. The potential at mj = 500 MeV is compatible 
with zero interaction. The analytic form ( |131| ) of the interaction potential 
was used as an input to the Schrodinger equation for a phase shift calculation 



according to scattering theory |36|. The resulting phase shifts are displayed 
in Figs. 26(b) and (d) for different values of the meson mass parameter used 
in the Schrodinger equation. 

One aim of this study was the comparison of results obtained from lattice 
QCD with results calculated from inverse scattering theory. There is a 



wealth of inversion results from hadron-hadron scattering |104]. We here 
restrict ourselves to cases with partial wave L = meson-meson scattering. 
The basis for any inversion is the existence of experimental data which, 
after a suitable interpolation and extrapolation, leads to a solution of the 
underlying ill-posed problem. For the tt-tt system the experimental situation 
is acceptable. The main source of experimental information for the inversion 
procedure is the analysis of the CERN-Munich experiment tt~p — > TT~TT~^n 



[]101| , 97, 105]. Even though there are ambiguities in the analysis of the vrvr 
final state, all phase shift analyses reach comparable results, which can be 
described by both meson exchange |^] and chiral perturbation theory (xPT) 
[p8| , 99 1 . This experimental situation is depicted in Fig. ^(f). The inversion 



results for several sources of phase shifts are shown in Fig. ^6Ke). While the 
isospin 1 = 2 channel is purely repulsive, the isospin 7 = channel shows 



an attractive component. A resonance is supported in this channel | 104 ]. 

The lattice mesons are coupled to isospin 1 = 2 but constructed without 
a specific angular momentum projection. Thus there is no one-to-one cor- 
respondence between the quantum numbers of the lattice mesons and the 



experimental data. The lattice potentials in Figs. 26(a) and (c) are attrac- 
tive with repulsion at short distances, whereas the potentials obtained from 
inverse scattering in Fig. ^(e) are only repulsive. This is expressed in the 
different form of the phase shifts in Fig. p6|(f) compared with those from 



the lattice QCD calculations in Figs. 26(b) and (d). One reason might be 



the incomplete projection to the correct quantum numbers. The construc- 
tion of good quantum numbers is not straightforward in the Kogut-Susskind 
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scheme. It is easier with the Wilson fermions but leads to extra terms in 
the correlator which have to be resolved numerically (see Appendix 0). 



To summarize Sees. 5.4 and we have performed two analyses partly 
in the same line both aiming at the extraction of a pion-pion potential from 
lattice QCD. The computation within staggered quarks suffered from the 
construction of spin observables whereas the more sophisticated calculation 
with Wilson fermions introduced a somewhat arbitrary cutoff at short dis- 
tances. As a general conclusion the reader has seen the current state of 
these trials which he might judge to be qualitative. Not only the increase 
of computer power but the efforts of physicists can make the results more 
accurate to compare with and predict experiments. 
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Figure 26: (a)-(d): Lattice results for meson- meson interaction potentials 
V{r) after the subtraction of the effects of the mirror particles, and phase 
shifts 5^ for quark masses "rnj = 100 MeV and mrif = 500 MeV. (e): vrvr 
L = 0, ^ 



/ = and 2 inversion potentials. 
[p8| , |99|] (dashed) and meson exchange § 
/ = and 2 phase shifts of Froggatt data 
solid line), Estabrooks et al. [100| (boxes). Grayer et al. jlOll] (diamonds). 
Manner |102|] (triangles) and Baillon et al. [103| (asterixes). 



Froggatt (sohd line), xPT 
(dotted), (f): from tttt L = 0, 
and interpolation (dots and 
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6 Conclusion and Outlook 



The field of hadronic physics has experienced a change of paradigm in recent 
years. From a modern perspective its primary goal now is the understanding 
of hadronic phenomena in terms of the fundamental degrees of freedom of 
the underlying theory, quantum chromodynamics |106|. Historically, phe- 
nomenological descriptions of hadron-hadron interactions have been with us 
for more than six decades while new questions related to the quark-gluon 
structure of hadrons have emerged as our insight has deepened over time. 
Thus hadronic physics now is closely intertwined with QCD. Dramatic ad- 
vances in computational technology, lattice field theory, and algorithms of 
computational physics, have moved unprecedented opportunities for funda- 
mental cognition within our reach. New accelerators and advanced detector 
design should of course also be mentioned in this context. Among the many 
topics tackled by lattice QCD - the hadron mass spectrum, hadronic struc- 
ture, heavy flavor and exotic hadrons, to name just a few - the physics of 
hadron-hadron interactions is one of the more challenging subjects. 

It is thus not surprising that the topic of this chapter is only in a nascent 
state. Concluding words with a quality of closure would be truly misleading. 
The current stage of hadron-hadron interactions on the lattice is driven by 
the tools of lattice technology with mass calculations being directly acces- 
sible from the euclidean formulation. All methods discussed in this chapter 
essentially use the energy spectra of certain systems as a starting point. 
Three ways to utilize those have been visited: finite volume methods and 
direct calculation of scattering phase shifts (Liischer's method), construction 
of an effective interaction from two-hadron states (perturbative potentials), 
and direct computation of the interaction in heavy-light systems (adiabatic 
potential). 

Of those methods the first one addresses measurable data directly, so 
far its most prominent application being aimed at scattering lengths. The 
other variants leave the realm of field theory in the hope of unveiling phys- 
ical mechanisms of the interaction. Perhaps, in the latter case the usual 
lattice technology constraints (related to chiral symmetry, quark masses, fi- 
nite volume, continuum limit, etc) are a lesser impediment to progress in 
the sense that qualitative insights into the physics of the interaction can still 
be gained with 'less than perfect' lattice parameters. 

Precision simulations reproducing experimental scattering data should 
probably not be the first priority at this time. To achieve this the field may 
not be advanced enough, particularly in the nucleon-nucleon case. Also, 
fermion methods dealing properly with chiral symmetry |107] should cer- 
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tainly be incorporated. 

Understanding hadronic interaction mechanisms, for example the impor- 
tance of gluon or quark degrees of freedom at certain relative distances, is 
probably more desirable at this stage. Emphasizing this aspect in future 
studies, for example heavy- light meson-baryon and baryon-baryon systems, 
seems a promising line of work that should shed some light on the QCD 
aspects of the hadron-hadron interaction. A related aspect is understanding 
of the origin of effective field theories [108]. 

Prom the viewpoint of lattice technology the computation of multi-quark 
correlation functions, which requires all-to-all fermion propagators, is cer- 
tainly the biggest problem. Considerable space was devoted to random 
source methods. Those are likely to remain pillars of hadron-hadron sys- 
tems simulations. In another development, not mentioned before, Baysian 
analysis methods for lattice generated time correlation functions promise 



access to excited states [109, |110| , an aspect crucial for the topic of this 
chapter. 

Finally, although it is true that the euclidean formulation of lattice QCD 
in imaginary time currently dominates the field, one should realize that it 
may not be the framework best suited to the task. It is not clear how to 
overcome the difficulty of describing asymptotic states, in the sense of LSZ 
scattering formalism, other than trying to solve the problem via spectrum 
computations at multiple lattice volumes 112 [. In the context of scattering 
the alternative hamiltonian formulation in real time is largely unexplored 

r2|. 
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7 Appendix A: Remarks on Staggered Fields 



The two most widely used schemes for discretizing fermion fields are com- 
monly known as 'Wilson fermions' |112] and 'Kogut-Susskind' or 'staggered 
fermions' 48, 114]. Wilson fermions are computationally more demand- 
ing than staggered fermions. While Wilson fermions explicitly break chiral 
symmetry, though they have a chiral limit which in practice can be recov- 
ered via a — > extrapolation of observables, staggered fermions retain 
a discrete remnant of chiral symmetry |115]. In the Wilson scheme con- 
struction of interpolating fields for hadrons with definite quantum numbers 
is transparent, and works more or less like in the continuum theory. This 
is opposite for staggered fermions, which entail complicated mathematics 
when constructing interpolating fields |115| , [73[ |. 

A case in point is an ordinary vr"*" interpolating field. First, the quark 
fields, Qf, are nonlocal operators which live on elementary hypercubes of the 



lattice |115, 114]- For a U(l) gauge model (no colors) we have 



qf{x) = V r^Uiv, x)xf{2x + r/) 



(A.l) 



built from the staggered fields Xf- The sum is over the set of the d- 
dimensional vectors r/ that point to the corners of the hypercube, r]^ = 0,1 
for fj, = 1 . . . d, with d = 4, see Ref. [114] for an odd number of dimensions. 
Further 



rr? = 7r---7d and U{7],x) = [Ui{2x)rK . .[Ud{2x + m ■ ■ -Vd-i)] 



Vd 



(A.2) 

and C, in (|A.l| ), is the dimension of the Clifford algebra. A vr^ interpolating 
field then is 

4>{x) = qd{x){j5 ® 75)9«(^) • (A-3) 

The direct product refers to row and column indices, respectively, of 
Clifford algebra matrices. Now, turning to the corresponding correlation 
function it can be shown that the latter reduces to a sum of local terms 

(^0^(f,t)</.(fo,to)) = (2^[2G(f,2t;xo,to)G*(x,2t;xo,to) 

+G{x, 2t + l;xo,to)G*{x, 2t + l;xo,to) 

+G{x, 2t - 1; xo, to)^* (x, 2t - 1; xq, to)]) . (A.4) 



For a concise exposition see Ref. |116| ]. This shows that the following equiv- 
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alent correlator may be used 

C(2)(t,to) = {Y,G{x,t;xo,to)G*{x,t;xo,to)) . 



(A.5) 



For our purposes a useful observation is that the local field 



(A.6) 



supplemented by a phase factor as shown gives rise to the correlator ( |A.5D . 
Hence, rather than using ( A.l ) and going through an involved derivation, 
we may as well work with the interpolating field 



+ ...+Xd^l+t 



Xd{x, t) Xu{x, t) 



(A.7) 



from the outset. The construction of correlators requires contractions be- 
tween the staggered Grassmann fields 



7t li . 

...Xf (x) Xf {x').^ 



6ffiG{x, x') 



(A. 



where n indicates the partners of the contraction, and G is the inverse 
fermion matrix, see (p6|). We also need 

... Xf (x) Xf (x )... = ...- Xf {x ) Xf (x) ... = ... - 5ffG{x , x) . . . . 

(A.9) 

In our application x' lives on a fixed time slice to whereas x runs over all 
lattice sites. This would require computing the entire inverse fermion matrix, 
as opposed to only one column. Here a symmetry can be employed. In the 
staggered fermion scheme the propagator satisfies |73] 

G{x',x) = {-lY' G\x,x'){-1)-'' with (-1)^ = (A.IO) 

which is actually obvious from (|2^). (The Wilson fermion version of this is 
Thus 



...If (x) Xf {x') ... = ...- 6ff{-ir'G\x, x'){-l)-- .... (A.ll) 
It is straightforward to check that the correlation matrices C^^^ and C^^^ of 



Sec. 3^ constructed with the fields ( [A.7| ) give exactly the same results as 
stated in ( p8| ) and ( p2[) re spectively. The reason is cancellation of the phase 
factors appearing in ( |A.7D and ( |A.ll ). 

For the purpose of concise presentation we omitted the phase factors in 
Sec. O. 
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8 Appendix B: Improved Lattice Actions 



The history of improved lattice actions originates with work by Symanzik 
1 117, 40 1 . The idea is to cancel discretization errors through adding addi- 



tional, carefully tuned, terms to the lattice action. Quantum effects (radia- 
tive corrections) play an important role in their construction. Symanzik's 
improvement program has received considerable attention |11^, |119|, |120|, 



121, |122| ] since its inception, but has only recently evolved into a form which 
is finding widespread application through the work of Lepage and Macken- 
zie [38|. Through a combination of beyond-elementary geometric structures 



in the lattice action and accounting for part of the quantum fiuctuations 
through numerically self-consistent renormalization lattice simulations now 
operate considerably closer to the continuum limit than ever before. With 
improved lattice actions results which can be compared to experiment are 
now within reach in many cases. With regard to the current topic (hadronic 
interaction) improved actions are important because they allow us to use 
large lattices, say L ~ 4 fm or so, with relatively moderate computing facil- 
ities. 

Improvement of the Lepage-Mackenzie fiavor has been well publicized 



[42|. Educational articles (for beginners) detailing the basics, ramifications. 



and physical significance of this development have been written by its origi- 



nators |123| , 124]. There is no point in duplicating those writings. However, 
given the importance of large lattices for hadronic interactions it is appro- 
priate to devote at least an appendix to the basics of this development. We 
will do so by way of selected examples. 

8.1 Appendix B.l: A Scalar Example 

In the lattice discretization of the space-time continuum derivatives are ap- 
proximated by finite differences. For example, considering a scalar field 
theory, 

Q2^^,^^Y,'^i^^^f^) + H^^-<^f^)-'^^(-) _ (B.l) 

(X 

In the classical limit, this approximation is 'good' if (p{x) is smooth on a scale 
larger than the lattice constant a. We may give a more precise meaning to 
this statement by saying that the Fourier transform 

^{p) = (27r)-2 J e^f-Xx) (B.2) 
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is zero (within some bound e) for momenta \p\ > vr/a. The error of the 
finite-difference approximation (BJ.) is o{a?)^ as a Taylor expansion shows. 
In fact we may systematically improve upon the classical discretization by 
including next-nearest-neighbor finite differences. Writing 

oo „ 
n=0 

we may solve the four equations 

g±a9^ = l±ad. + -a^dl±-a^dl + —a^df,±... (B.4) 
^ 2 ^ 6 ^ 24 ^ 

^±2aa, ^ l±2adf, + 2a^dl±U^df, + la'^df,±... (B.5) 
for d^, n = 1 . . . 4, m terms of powers of single-step shift operators 

Specifically we find 

ad, = ^(AW-A(-))-l(A(+)2-A(-)2) + o(a5) (B.7) 
a'dl = l(AW + A(-)-2)-l(AW2 + A(-)2-2)+o(a^).(B.8) 

Introducing the finite-difference derivatives 

a{+) a (-) a(+) _l a(-) 9 

A« = ^^^— ^ and A(^) = ^^ + (B.9) 

we may also write 

d, = AW-la2(AWA(2)+A«A(2)) + o(a^) (B.IO) 

dl = ^A(2)-^A«2 + o(a^). (B.ll) 

Note that A^^'* and A^^'* commute. Also, it is straightforward to generalize 
the expansions to higher orders in a in a systematic way. Thus the continuum 
action, say for a scalar Klein-Gordon field with quartic self interaction 

S[<P] = J d'x[-^{d,<pf + imV' + ^^P'] (B.12) 
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may then be discretized, alternatively, as 



(B.13) 



or 



2 ^ 



3 ^' 3 y ; 



(B.14) 



The above lattice actions are equivalent up to o(a^). Evidently the advan- 
tage of o(a^) improvement is that one works close to the continuum limit 
(a — > 0), or from a numerical point of view, at the same 'closeness' to the 
continuum theory on can afford a coarser (larger a) lattice which translates 
into less computational cost. 

Since a classical field 4>{x) is smooth on some scale, say a, improvement 
should be expected to work at some level of accuracy. 

However, upon quantizing the theory, smoothness of the field is lost. In 
the generating functional 



Z[J] 



d(/)] e 



(B.15) 



(j){x) and 4>{y) are independent integration variables for x ^ y even if x 
and y are infinitesimally close. At a given scale, say a, those fluctuations 
\4>{x) — (p{y)\ can in principle be large for |x — ?/| <C a. Their importance 
is determined by the Boltzmann weight in the euclidean path integral. In a 
manner of speaking we say that the quantum fields are rough on any scale, 



see Fig. 

This means that the Fourier transform ( [B.2D is unbounded and, generally 
speaking, momentum-space integrals diverge. The well-known solution is 
renormalization. In a nutshell, introducing a momentum cutoff A, which 
makes the momentum-space integrals finite, then renormalizing the fields 



4>{x) = Z^^(I)a{x) 



(B.16) 



couplings, etc, and, at the level of physical quantities (observables) removing 
the cutoff A, becomes an integral part of the quantum field theory consid- 
ered. 
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Figure 27: Roughness of the quantum field on any scale. The smooth curve 
illustrates a Fourier average with a finite momentum cutoff A. 



Clearly, lattice discretization introduces a natural cutoff momentum 

A = -. (B.17) 

a 

In a way, the roughness of the fields at scales less than the lattice constant a is 
rendered harmless at the expense of introducing renormalization. (One may 
think of Z\, for example, as playing a role similar to a dielectric constant 
for a medium of electric dipoles. The small-scale fluctuations of the electric 
field on the scale of the dipoles are lumped into Ke-) 

The field renormalization constant is without consequence in the present 
example because it can be made to disappear by changing integration vari- 
ables 4>a{x) in the generating functional (|B.15| ). Thus 

Z[Ja] = const j [d</>A] e-^A[<^A]+JA-<AA^ (b.18) 



where, S\ is the same as ( B.13 , B.14| ) except that /3a = PZ^'^ and 



gZ~j^ . Improvement does not change the physics of the quantum theory in 
this particular example of the scalar interacting Klein-Gordon field. 

This culture of thinking has, perhaps, impeded a much earlier discovery 
of the type of improvement now routinely used in lattice gauge simulations. 
The above line of thought does lead to a genuinely different situation in the 
case of gauge theories. 
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8.2 Appendix B.2: Improvement of a Pure Gauge Theory 



We start from the Schwinger representation 125 | of the (lattice) hnk vari- 
ables 



U,i{x) = V exp (^g j A^{x + iji) 



(B.19) 



where V denotes path ordering. The latter ensures proper behavior under 
gauge transformations. The continuum limit of the classical, smooth, field 
is obtained by expanding U^{x) into a power series of a. Some formalism 
proves useful. Motivated by 



/ diA^{x + 5,fi)= [ dCe^^^A^ix) = aTi{ad^)A^{a 
Jo Jo 

we define the set of (integrated) translators 

oo 

Tn{X) = J2 



k=0 



{k + n)\ 



Those satisfy 

To{X) = 



Tn{X) = 1 



1 



n + 1 



Tn+l{X)X. 



(B.20) 



(B.21) 



(B.22) 



A straight path along m links in direction fi involves 

/•ma 

/ dCAuix + Cfi) = maTi{mad^)Au{x). (B.23) 
Jo 

These can be used as building blocks for various loop operators, for example, 



the planar m x n plaquettes, see Fig. 28, 



-ReTrT'exp 
3 



-ig 



pma 

Ig d^Af,{x + ^fL) 
Jo 

pna 

/ A,y{x + mafi + ^u) 
Jo 

/•ma 

/ d^ A^ {x + naO + ^/i) 

^0 

/•na 

/ d^ Ay{x + iv) 
Jo 



-ReTrT^exp \^gmna?'T\{nadv)T\{madp) 
{d^A,{x) - d,A^{x))\ . 



(B.24) 



(B.25) 
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A 

x+a V 




Figure 28: Geometry of some loop operators, elementary 4-link plaquette, 
6-link planar rectangles, bent rectangle, and parallelogram, considered for 
improved gauge field actions. 



The combination of path ordering V and taking the color trace Tr ensures 
that P^™"""^ is gauge invariant. Manifest gauge invariance may be recovered 
by replacing with the covariant derivative 

D^ = d^- igA^{x) . (B.26) 

As a consequence, the expansion of ( B.25| ) into a power series of a only con- 
tains gauge- invariant terms. Specifically, for the 4-link elementary plaquette 
we obtain 

-^a^Tr {gF^,{x){Dl + Dl)gF^,{x)) + (.(a^) (B.27) 
and for the 6-link planar rectangle 



2 
— ( 
3 

^ a^Tr {gF^,[x){ADl + Dl)gF^,{x)) + 0(0^) (B.^ 



18 



72 



For P^J,""^^ we have to replace ADl + Dl with Dl + AD^ in (|B.28D . Thus the 



classically improved Wilson action for the pure gauge theory is 



x2)^ 



(B.29) 

with C2 chosen such that the terms of the classical limit cancel. This is 
achieved by 

C2 = — — (classical) 



in which case f3 = and 
Sg[U] ^ 



X ii<u 



/?«^EE 

X f. 

E 



10 



Tr(F^,(x))2+o(a3) 



i^Tr(F^,(x))2+o(a3 



(B.30) 

(B.31) 
(B.32) 



In the quantized theory expansion of the gauge fields into a power series may 
no longer be justified. Nevertheless, we expect the expansion to be 'good' 
if the expectation values of the plaquette operators (B.27) and ( B.28| ) are 
'close' to one, 

X ii<u 



This condition can be easily tested in practical situations. In Fig. 29 for 
example, are shown the lattice averages of the 4-link plaquette and the 
planar 6-link rectangle. The deviations from one are sizeable and increase if 
the gauge field coupling {g'^ oc f3~^) is made larger. This effect is, of course, 
due to quantum fluctuations. Its significance is that 

• Quantum fluctuations spoil the classical continuum limit. 

Again, the standard procedure to deal with short-ranged (ultraviolet) fluc- 
tuations is renormalization. An ultraviolet momentum cutoff is in place 
through lattice discretization. The gauge fields are represented by the link 
variables ( B.19 ), thus replace 

U^{x) ^ u-'U^ix) (B.34) 

where u is a renormalization constant. Note that without the 6-link terms 
in the action S'G'[f7], see ( B.29 ), the factor u would just renormalize the 
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0.0 ^ 

4 5 6 7 8 9 10 11 12 13 



Figure 29: Expectation values of the 6-link plaquette (planar rectangle), the 
4-link (elementary) plaquette, and the corresponding tadpole renormaliza- 
tion factor u for versus /?. 

gauge field coupling, and be harmless (and useless) very much like in the 
example of Sec. ^.11 . Link variable renormalization does not do anything 
if only the standard Wilson action is used. The crucial point here is that 
the 6-link improvement terms in the action carry two more factors of link 
variables. Their fluctuations (in the mean- field sense) can be divided out by 
two renormalization factors. Hence, using 

= ~ (B.35) 

in the action results in the cancellation of the o(a^) terms, as above, between 
the 4-link and 6-link terms in the action. 

• Classical improvement and gauge link renormalization conspire to can- 
cel discretization errors in the quantized lattice field theory. 

This state of affairs is characteristic for a gauge theory. In the lingo of 
perturbation theory the classical limit is known as the tree-level approxima- 
tion, because no loop diagrams are taken into account. The renormalization 
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factor u comprises the effect of so-called tadpole diagrams, which are singly- 
attached loops to the gluon propagator. 

To be sure, there are radiative corrections to the gluon propagator which 
are not contained in the tadpole diagrams. Some of those induce o(a^) 
corrections oc agin /a) and, containing loops, are harder to come by [122| , p] ]. 

The tadpole factor, on the other hand, can be monitored during an actual 
simulation. Starting with an arbitrary value, say u = 1, after a number of 
iterations we simple 'measure' the elementary plaquette (for example) and 
redefine 

- = (^EE^£"V/^ (B.36) 

The new value is then used in the action SgI^^] and the iterations continue. 
Another u is then 'measured' with the new action, and so on. In this way 
the 

• tadpole renormalization factor is determined self consistently during 
the simulation. 

Numerical work has shown that tree-level and tadpole improvement 
(even without loop corrections) already account for much of the discretiza- 
tion errors The effect of improvement is quite dramatic. A well- 
known example [124] is the restoration of rotational invariance of the static 



quark-antiquark potential, usually parametrized as V{r) = Kr — 7r/12r + C. 
On- axis distances r /a = 1,2, 3, 4... and off- axis distances r/a = \/2, \/3, \/5 . , 
which are supposed to fall on a smooth curve V{r), are off by as much as 
38%. With a tree- level and tadpole improved action errors are reduced to 
the one-percent level! 

Improvement schemes are not unique. Various loop geometries (bent, 
twisted rectangles, etc) may be used to cancel discretization errors. A cor- 
responding statement of course holds for improving the fermionic part of the 
action. 

8.3 Appendix B.3: Improvement of a Fermionic Action 

A simple fermion action, which is improved up to the same order as the 
gauge field action, can be constructed by employing next-nearest-neighbor 
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terms \12t, 44, |127| ]. Consider the following 



S 



'(0) 
F 

'(1) 



^(2) 



a^^'tp{x)ip{x) (B.37) 

X fl 

+Tp{x + ail){r + j^)Ul{x)iP{x)\ (B.38) 
a^J2Yl - lt^)U,^{x)Uf,{x + afi)ij{x + 2afi) 

X fl 

+i^{x + 2afL){r' + -ff,)Ul{x + afi)Ul{x)ip{x)\ (B.39) 

where r, r' are parameters. Knowing about the essential ingredients for tree- 
level and tadpole improvement we construct a fermionic action as a linear 
combination 

(B.40) 



Sf — coSp^ + —Sp^ + -^Sp^ . 



u " W 

The factor u is the tadpole renormalization ( [B.36D which reflects the relative 
number of gauge links in Sp^ and Sp^ , Sp^ . We need the following classical 



continuum limits 



S 



r{2) 



8r a^^ipip- ^ ipliiD^ip (B.41) 

X X 
X X 

8ra^ ^tpip- Aa'^ ^ i'lfiDf^il^ (B.42) 

X X 



Remarkably it is possible to simultaneously cancel the o(a^) and o(a^) terms. 
Cancellations must be realized on the tree level, u = 1. The requirement is 





A nontrivial solution exists if 

r' = 2r with 



ci + 8c2 = . 



(B.43) 



(B.44) 
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6 



We solve the last equation writing 

4 , 

Cl = --C C2 

where c' is a normalization. Still on the tree level, we then have 
Sf {co-8c'r)a^^ip7p + 2c'a'^^ipj^D^iP + o{a 

X X 

Co — 8c'r 



(B.45) 



2da 



(B.46) 



for a — > 0. Thus, in the continuum limit, the bare fermion mass is = 
(co — 8c'r)/2c'a. The choice 

r = 1 (B.47) 

for the Wilson parameter takes care of the doubling problem (violating chiral 
symmetry). In order to meet with the usual conventions we introduce the 
'hopping' parameter k through 



and choose the normalization 



c = KUCq 



Co = ^ 



(B.48) 



(B.49) 



Thus, finally, the improved lattice action ( B.40 ) becomes 



X fj, 



X fl 



+xi;{x + 2aA)(2 + -/f,)Ul{x + a/i)[/t(x)V'(x) 



(B.50) 



With this improved action we are hopefully 'close' to the continuum limit, 
so the fermion mass is 'close' to mj = (1 — 8Ku)/2Kua with a critical value 
of He = l/8n for = 0. 
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8.4 Appendix B.4: More on Highly Improved Actions 



An example of a dispersion function computed with the o{a^) improved 
action ( |R29| , p35| ) and ( p3o| ) is shown in Fig. |0[ The fihed circles are the 
energies of a pseudoscalar meson with momentum p from a 5^ x 7 lattice 
at P = 5.0. The data are close to a smooth (fitted) curve, as opposed to 
the open circles, which are from the standard lattice dispersion relation for 
bosons 



12 



E^=M^+P^-PV36 




iioL improved 



1//c=0.174 M = 1.42 



15 



Figure 30: Energy-momentum relation for a pseudoscalar meson obtained 
from a lattice simulation with the improved action (filled circles), and a 
fitted curve. Also shown are the unimproved lattice dispersion relation for 
a free Bose field IS^l, and the continuum relation. 



Another example of 'closeness to the continuum' achievable with im- 



proved lattice actions is shown in Fig. |31|. The nucleon-rho mass ratio is 
plotted versus the rho mass. Open plot symbols correspond to the standard, 
unimproved, Wilson action (gauge and fermion fields). The filled symbols 
refer to various 'flavors' of improved actions. The experimental value of 
mN/mp is 1.22. The two curves in Fig. meant to 'guide the eye', illus- 
trate dramatically the 'quantum leap' provided to us by improved actions: 
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Physical results, close to the continuum limit, can be obtained from simula- 
tions on relatively coarse lattices with lattice constants a about a factor of 
~ 4 larger compared to the standard Wilson action. Systematic errors still 
are large though, at this time. 

2,2 I , , , ^ , , 1 




1.0 ' ' ' ' ' ' 1 

.0 .1 .2 ,3 A 

a [fill] 



Figure 31: Nucleon p-meson mass ratio versus the lattice constant a. The 



scale for a was obtained from Mp = 770 MeV. See Ref. |28|, |T2|, |T3g] to 
trace the references for the data points. 



The lattice actions discussed in Sees. and |8^ are improved up to two 
leading orders in the lattice constant. Those actions are commonly referred 
to as o(a^) improved, or highly improved actions. It should be mentioned 
that highly improved actions are not without problems. The dispersion 
function E{p) for free massless quarks, as calculated from the zeros of the 
inverse lattice propagator, exhibit unphysical branches that describe addi- 
tional particles, or ghosts. Hopefully their masses are sufficiently large, so as 
not to interfere with the propagation of the physical quarks on the lattice. 
The effect of those particles reveals itself in the early-time-slice behaviour of 
hadronic correlation functions. Thus, in practice, their effect can be isolated 
and hardly poses a problem. 

A class of highly improved actions that has received much attention 
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recently is known as D234 actions, the nomenclature taken from the way 
derivative operators are used |131]. Here, in addition to the techniques dis- 
cussed above, one utilizes a construct known as the Sheikholeslami-Wohlert 
(SW), or clover, action |43] as well as anisotropic lattices. The latter have 
different lattice constants and at for space and time respectively, typically 
as > at- The purpose of combining those features is to push doublers, or 
ghosts, into the high- mass region where they can do little harm to hadronic 
physics. Actions of this class are well documented, technical issues are in 
the foreground at this time [131, 132, 133, 134 1. We refrain from discussing 
these further. 



8.5 Appendix B.5: Clover Leaf Fermion Action 

The SW action deserves being mentioned here because of its widespread use 
and relatively simple implementation. 

The idea is to consider (integration) variable transformations of the 
fermion fields in the path integral that defines the lattice QFT. This will 
of course not alter the physics of the theory but, in general, change the 
appearance of the lattice action written in terms of the transformed fields. 
The additional terms are called redundant. Sheikholeslami and Wohlert [^] 
considered (among an extensive list of other possibilities) the transformation 

vn — — TO 

i; = [l- -(7mA^ - m'M ^ = V^'[l - -(A,,7^ - m')] (B.51) 

where ip\^' refer to the original field variables and mass m' . Working out 
the Dirac action we obtain 

— — 7*Qj 1 

V;'(7^ + mOV'' = V'[7mA/. + - Y (A/. A^ + -(J^^gF^^)]ilj + o{a^) (B.52) 

with ^ 

m = m'(l + —ram!) . (B.53) 

The a ■ F term stems from the identity 

{Df,jf,f = Df,D^ + ^af,ugF^u with cr^^ = ^[7^,7,.] . (B.54) 
This leads to the SW, or clover, action 

X fl,U 
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Figure 32: Geometry of the paths used in P^iy{x) for the clover action. 

which is an additive part to the standard Wilson fermion action. The name 
clover action is apparent from Fig. ^ which shows the geometry of link 
variables involved in 

P^^(x) = ^(^ Uf,{x)U^{x + ail)Ul{x + aiy)Ul{x) 

— Ul{x — av)Uj^{x — afi — av)Up{x — afi — ai))U^{x — ajl) 
+ Ui,{x)Uj^{x — afi + ai))Ul{x — ajl)U^{x — afi) 

—U^{x)Ul{x + afi — ai>)Ujj^{x — ai')Uu{x — av)^ . (B.56) 

The action is o(a) improved on the tree level, and local, meaning that no 
'hopping' to next neighbors or beyond occurs. This feature is computation- 
ally desirable since it leads to simpler and faster codes. 



The clover action and the action discussed in Sec. 8.2 are related [135]. 
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